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Abstract 

An effective field theory approach is used to describe quantum matter at greater-than- 
atomic but less-than-nuclear densities which are encountered in white dwarf stars. We focus 
on the density and temperature regime for which charged spin-0 nuclei form an interacting 
charged Bose-Einstein condensate, while the neutralizing electrons form a degenerate fermi 
gas. After a brief introductory review, we summarize distinctive properties of the charged 
condensate, such as a mass gap in the bosonic sector as well as gapless fermionic excitations. 
Charged impurities placed in the condensate are screened with great efficiency, greater than 
in an equivalent uncondensed plasma. We discuss a generalization of the Friedel potential 
which takes into account bosonic collective excitations in addition to the fermionic excita- 
tions. We argue that the charged condensate could exist in helium-core white dwarf stars 
and discuss the evolution of these dwarfs. Condensation would lead to a significantly faster 
rate of cooling than that of carbon- or oxygen-core dwarfs with crystallized cores. This pre- 
diction can be tested observationally: signatures of charged condensation may have already 
been seen in the recently discovered sequence of helium-core dwarfs in the nearby globular 
cluster NGC 6397. Sufficiently strong magnetic fields can penetrate the condensate within 
Abrikosov-like vortices. We find approximate analytic vortex solutions and calculate the 
values of the lower and upper critical magnetic fields at which vortices are formed and de- 
stroyed respectively. The lower critical field is within the range of fields observed in white 
dwarfs, but tends toward the higher end of this interval. This suggests that for a significant 
fraction of helium-core dwarfs, magnetic fields are entirely expelled within the core. 
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Introduction and Summary 



It is an everyday experience that by simply changing the temperature of a substance we 
can abruptly and dramatically change its macroscopic properties. At normal earthly tem- 
peratures and pressures we observe, e.g., solids, liquids and gases, and the phase transitions 
between these states. As we consider more extreme conditions out of the realm of our daily 
experience, say at temperatures near absolute zero or at the high densities that exist in the 
cores of many stars, phase transitions continue to occur and new states of matter emerge. 
The principles of fundamental physics enable us to predict the properties of states of matter 
that have yet to be observed. 

At low temperatures and high densities, quantum mechanics becomes essential for de- 
scribing the properties of a state of matter. An example of this is the Bose-Einstein conden- 
sate. Below a certain critical temperature the thermal de Broglie wavelengths of an ideal (or 
nearly ideal) gas of bosons will begin to overlap. This critical temperature corresponds to 
when the quantum-mechanical uncertainties in the positions of the particles becomes greater 
than their inter-particle separation. Below this temperature, the statistics governing a gas 
of indistinguishable bosons dictate that the particles will "condense" into the same quantum 
state. The first gaseous Bose-Einstein condensate was created in a laboratory seventy years 
after its existence was first predicted [1]. However, the extreme conditions required for the 
existence of such a quantum substance can occur outside the lab as well, in astrophysical 
objects. 

The cores of white dwarf stars are composed of a particularly dense system of nuclei and 
electrons, with an average inter-particle separation much larger than the nuclear scale but 
much smaller than the atomic scale. Because white dwarfs have exhausted their thermonu- 
clear fuel, they evolve by cooling. At high temperatures, the equilibrium state of the system 
is a plasma. As the system cools below a certain temperature, the energy of the Coulomb 
interactions will significantly exceed the classical thermal energy. Then, the standard clas- 
sical theory holds that the nuclei will arrange themselves in such a way as to minimize their 
Coulomb energy, forming a crystal lattice [2]. It is expected that in most white dwarfs, 
consisting of carbon, oxygen, or heavier elements, the crystallization transition takes place 
in the process of cooling (for a review, see, e.g., Rcf. [3]). 

In this work we will argue that for white dwarf stars with helium cores quantum effects 
become significant before the crystallization temperature is reached. In these cases, the de 
Broglie wavelengths of the nuclei begin to overlap before crystallization can occur. Then, 
because the helium-4 nuclei are bosons, the quantum-mechanical probabilistic "attraction" 
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forces the nuclei to undergo condensation into a zero-momentum macroscopic state of large 
occupation number. The nuclei minimize their kinetic energy, while their collective fluctua- 
tions (phonons) have a mass gap. The majority of the phonons cannot be thermally excited 
since the phonon gap ends up being greater than the corresponding temperature. There- 
fore, after the phase transition all the thermal energy is stored in the near-the-fermi-surface 
gapless excitations of quasi- fermions. We refer to this state as a charged condensate [4-9] . 

Because the repulsive Coulomb interactions between the ions dominate over the thermal 
energy of the system, the bosonic sector of the charged condensate is strongly coupled. This 
is in contrast to the usual neutral Bose-Einstein condensate. In this work we summarize an 
effective field theory approach to describing the charged condensate. Within this framework 
we find properties of the charged condensate that are distinct from its neutral counterpart. In 
particular, as mentioned above, we find that the spectrum of the collective bosonic excitations 
is gapped and the bosonic contribution to the specific heat is exponentially suppressed at low 
temperatures. Then, most of the entropy of the system is stored in the near-the-fermi-surface 
gapless fermionic excitations. 

Furthermore, we find that electrically charged impurities in the condensate are screened 
to a high efficiency, more effectively than in an equivalent uncondensed plasma. The static 
potential contains an exponentially suppressed term as well as a long-range oscillating piece. 
The latter is due to gapless fermion excitations, and is similar to the Friedel potential. 
However, the potential is also suppressed due to an attractive phonon interaction, and we 
obtain an expression which has a long-range oscillatory nature but is highly suppressed 
compared to the conventional Friedel potential. 

Such properties of the charged condensate have consequences for helium white dwarfs. 
A condensed core dramatically affects the cooling history of the helium white dwarfs - they 
cool faster than those with crystallized cores. As a result, the luminosity function exhibits a 
sharp drop-off below the condensation temperature [9] . Such a termination in the luminosity 
function may have already been observed in a sequence of the 24 helium-core white dwarf 
candidates found in the nearby globular cluster NGC 6397 [10]. 

While the focus of this paper is the application of the charged condensate to helium 
white dwarfs, our methods are general and can be readily apphed to other systems. The 
bosonic field can be generalized to any fundamental scalar field or a composite state, and 
the electromagnetic interaction could be replaced by any U{1) abelian interaction. The 
applications of a new state of matter are potentially diverse. 

The structure of this work is as follows. In chapter 1 we review the condensation of a 
gas of neutral bosons and give the expression for the standard critical temperature at which 
condensation occurs. We review the condensation of a gas of weakly interacting bosons in 
the formalism of both a relativistic and a non-relativistic effective field theory. 

In chapter 2 we give arguments for the existence of the charged condensate. A general 
mechanism for charged condensation in the context of a relativistic field theory is presented 
there and the spectrum of small perturbations above the condensate is determined. Fur- 
thermore, it is shown in chapter 2 that electrically charged impurities are screened to a high 
degree due to an attractive phonon interaction. There we consider the effects of fermion exci- 
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tations on the electric potential and derive a generalized Priedel potential for the condensate. 
We also briefly discuss the generalization of the Kohn-Luttinger potential. 

In chapter 3 we argue that in helium-core white dwarf stars, the helium-4 nuclei may 
condense as they cool, instead of crystallizing. A low-energy effective field theory description 
of the helium-4 charged condensate is developed, from which we recover the same charac- 
teristic properties that we found in the relativistic theory. Furthermore, we consider the 
cooling rate of dwarf stars with condensed cores and show that the age of such dwarfs would 
be significantly shorter than those with crystallized cores. 

In chapter 4 we look at the magnetic properties of the charged condensate. One would 
expect these to be similar to type II superconductors. Indeed, we find vortex-type solutions 
for magnetic flux tubes in the charged condensate. Prom these we determine the magnitude 
of the external magnetic fleld for which it becomes energetically favorable to form vortices. 
We discuss the applicability of the vortex solutions to magnetized helium-core white dwarfs, 
and also consider the effect of a constant rotation on the magnetic fleld in the condensate of 
helium-4 nuclei. 

Let us make a few comments on the literature to emphasize the differences of the present 
approach. The condensation of non-relativistic charged scalars has a long history, the original 
works being those by Schafroth [11] in the context of superconductivity, and by Foldy [12] in 
a more general setup. An almost-ideal Bose gas approximation was assumed in those studies. 
For this assumption to be valid, densities had to be taken high enough to make the average 
inter-particle separation shorter than the Bohr radius of a would-be boson-antiboson bound 
state [12]. If the fermion number density is denoted by Jo, and the mass of the scalar by ttih, 

1 /3 

this would be the case if Jg ^ ctem'mH- However, for a helium-electron system the above 
condition would translate into super-high densities, at which nuclear interactions would 
become signiflcant. Instead, in this work we study charged condensation in the opposite 

1 /3 

regime, Jq <^ acmfnH, where the nuclear forces play no role. As a result, certain properties 
of the system - such as important details of the spectrum and the screening of electric charge 
- are different. Also, our method, which is based on symmetry and fleld theory principles, 
is different. 

In the context of a relativistic fleld theory the condensation of scalars was discussed in, 
e.g., Refs. [13-15]. Pion condensation due to strong interactions is well known [13]. In 
our work strong interactions play no role (except for providing the nuclei). It was shown 
in Ref. [14] that a constant background charge density strengthens spontaneous symmetry 
breaking when the symmetry is already broken by the usual Higgs-like nonlinear potential 
for the scalar. In our work the scalar has a conventional positive-sign mass term and no 
Higgs-like potential. The fact that the conventional-mass scalar could condense in a charged 
background was flrst shown in [15], in the case that the scalars have a nonzero chemical 
potential (see brief comments after eq. (4.6) in [15]). Moreover, a somewhat similar bosonic 
spectrum was already discussed in a work [16] in the context of superconductivity. However, 
the dominant contributions to the thermodynamics of the charge condensate discussed here 
are due to near-the-fermi-surface electron excitations, which were absent in the the system 
considered in Refs. [15,16]. 
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The possibility of having a charged condensate in hehum white dwarfs was previously 
pointed out in Ref. [17], where the condensation was studied using an approximate variational 
quantum-mechanical calculation in conjunction with numerical insights in a strongly-coupled 
regime of electromagnetic interactions. The degree of reliability of such a scheme is hard to 
assess. Furthermore, using the ordinary neutral Bose-Einstein (BE) condensation to describe 
the charged condensate, as is done in a number of works in the literature, is hard to justify. 
As we will see, properties of a neutral BE condensate differ significantly from those of the 
charged condensate considered here. For instance, the specific heat at moderate temperatures 
in the former is due to a phonon gas, while in the latter it is due to the degenerate electrons. 
We will discuss this feature of charged condensation and others in greater detail in what 
follows. 

The novel feature of our work is the development of an effective field theory approach 
to condensed matter at greater-than- atomic and less-than-nuclear densities, and its appli- 
cation to helium-core white dwarfs. This field-theoretic framework allows us to study the 
more subtle aspects of the condensate, including its spectrum and properties of its magnetic 
vortices. The latter is especially hard to analyze without the field theory approach. The 
generahzation of the Friedel potential and the Kohn-Luttinger effect to a system with col- 
lective excitations of both bosonic and fermionic nature was obtained, using this method, in 
our work [7]. In addition, we are able to make concrete predictions as to the fast cooling of 
helium-core white dwarfs. 

In Ref. [18,19], A. Dolgov, A. Lepidi, and G. Piccinelli have performed a one-loop calcu- 
lation and found finite temperature effects in a general setup with condensed bosons. These 
authors also obtained infrared modifications of the static potential for the scalar 
we had earlier. Our approach and that of Refs. [18, 19] are complementary. We empha- 
size understanding the charged condensate in terms of effective field theory and low-energy 
collective excitations. 

Besides the first chapter which introduces the field theory description of a neutral con- 
densate, the bulk of the present paper is based on our previous works on charged conden- 
sation [4-9]. However, it also contains new technical and conceptual details that have not 
been published elsewhere. 



Notations and Conventions: We work in natural units where H = c = = 1 unless 
explicitly stated otherwise. The signature of the metric tensor taken to be — , — , — ). We 
use Heaviside-Lorentz units for Maxwell's equations. Accordingly, the fine-structure constant 
is given by 

1 

(y — — ~ 

47r " 137 ■ 
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Chapter 1 

The Neutral Condensate 



1.1 statistical mechanics of condensation 

At low temperatures and high densities, a gas of bosons exhibits fundamentally different be- 
havior from a gas of fermions. When the concentration of particles is sufficiently high so that 
the thermal de Broglie wavelengths of the particles begin to overlap, then quantum effects 
become important; the particles must be treated as truly indistinguishable with fermions 
obeying Fermi-Dirac statistics and bosons obeying Bose-Einstein statistics. Because bosons 
do not obey the Pauli exclusion principle, they are free to occupy any state of the system 
in arbitrarily large numbers. Thus at sufficiently low temperatures, the ground state of 
a system of bosons will be macroscopically occupied, forming a Bose-Einstein condensate. 
In what follows we briefly review the statistics of a gas of bosons that lead to the critical 
temperature Tc at which Bose-Einstein condensation occurs. 

For a gas of bosons, the average number of particles in a state of energy e is given by the 
usual Bose-Einstein distribution: 

NMe,T)^ ^^^_^y^_^ . (1.1) 

Let us consider a non-relativistic gas of free bosons with energy e = p^/2m, where m is the 
mass of the particle. Take the ground state of the system to have zero energy: eo = 0. Then, 
as the non-relativistic chemical potential^ approaches zero fi^R 0~, it is clear from the 
above expression that the number of particles in the ground state will become large: 

T 

No as iiNR 0" . (1.2) 

The total number of particles in the system is the sum of the number of particles in the 
ground state and in the excited states: Nff^t — Nq + N^^- In the ^nr limit, one can 
calculate the number of particles that remain in excited states. It is just the integral over 

^The non-relativistic chemical potential is related to the relativistic chemical potential by ^inr = /U — m. 
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all momentum states of excited particles: 



Here V is the volume of the system. From this expression we can define a critical tem- 
perature Tc at which the number of particles in excited states is equal to the total number 
of particles. Above this temperature, there are negligibly few particles in the ground state 
and a nonzero (negative) chemical potential must be restored in expresion (1.3). Below this 
critical temperature the ground state will be macroscopically occupied. Taking N^x = -/Vtot 
and the density of particles to be n = Nfot/V we find the usual critical temperature: 

2./j^y/» (14) 

Here ( is the Riemann zeta function: ^(|) ?a 2.612. If we take the interparticle separation 
to be d = (|7rn)~-^/^ the critical temperature becomes 

1 27 

^ ^ . (1.5) 

This critical temperature corresponds to when the thermal dc Broglic wavelengths of the free 
particles overlap with each other. In other words, the condensation begins to occur when 
the quantum mechanical uncertainties in the positions of the particles become greater than 
the interparticle separation. 

The critical temperature is remarkable in that it can be significantly higher than the 
energy of the first excited state of the system. Classically, we would expect that, in a 
system at temperature T, most particles would be in single-particle states with energy of 
order T, for arbitrarily small T. However, the quantum statistics of a gas of bosons at 
low temperatures lead us to a strikingly different conclusion. At low temperatures, but 
at temperatures still much higher than the energies of the lowest accessible excited states, 
bosons prefer to macroscopically occupy the ground state. 



1.2 The effective field theory 

For a large number of quanta that form a macroscopic state, the state can be adequately 
described in terms of an effective field theory of the order parameter, and its long wavelength 
fiuctuations. In particular, we will look for classical solutions of the equations of motion of the 
effective order-parameter Lagrangian. How can a classical solution describe the condensate 
which is an inherently quantum phenomenon? Denote the particle creation and annihilation 
operators by and Oq respectively; then the quantum-mechanical noncomutativity of these 
operators, alaQ—aoal ~ h, becomes an insignificant effect of order 0{h/N), when the number 
of particles in the condensate state, {alao) ~ N, is large enough, ^ 1. Thus, the classical 
description of a coherent state with a large occupation number - the condensate - should be 
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valid to a good accuracy [21]. On the other hand, collective excitations of the condensate 
itself should be quantized in a conventional manner. 

The above arguments lead to the following decomposition of the order-parameter operator 
describing the condensate: 

$ = $^ + 5$, (1.6) 

where $cz denotes just a classical solution of the corresponding equations of motion, and 
describes the condensate of many zero-momentum particles, while 5^ should describe their 
collective fluctuations. 

In what follows we focus on the zero-temperature limit, even though realistic tempera- 
tures in, say, helium white dwarfs arc well above zero (for calculations of the finite tempera- 
ture effects, see [18, 19]). We will justify the validity of the zero-temperature approximation 
as we proceed (see, in particular, section 3.2). 

1.2.1 The relativist ic EFT 

To describe the neutral condensate of bosons in terms of an effective field theory, we adopt 
the simplest Lagrangian for the order parameter that exhibits the condensation. The field 
is a complex scalar field with a right sign mass term > and a repulsive self-interaction 
with interaction strength A: 

jr=\d,<P\'-mU*<l)-\{<j>*^f. (1.7) 

This Lagrangian could contain higher order terms, however, they are generally suppressed 
by the short-distance cut-off of the theory and are irrelevant for our considerations. Such a 
relativistic model of condensation was considered in Ref. [22]. The first microscopic theory 
of condensation in weakly interacting bose gases was developed by Bogohubov [23] . 

For convenience, we switch notation and write (f) in terms of its modulus and phase: 
(f) — :^cre*°. The Lagrangian (1.7) becomes: 

^ - lid.ar + lid^afa^ - \mla^ - . (1.8) 

Written in this form, it is evident that a nonzero value for d^a acts as a tachyonic mass for 
the scalar.^ 

Varying the Lagrangian w.r.t a gives the conservation of the scalar current density: 

d^[{d^a)a^] =0. (1.9) 

If we take djO. = 0, then the number density of scalars is constant in time. Let us denote 
this number density by Jq = dgaa^. 

^Onc could also introduce a chemical potential for the scalar which would have the same effect, however 
this term can be effectively absorbed into a redefinition of doa. 
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Varying the Lagrangian with respect to a gives the following equation of motion: 

Ua=\{d^cif -m\\a -\a^ . (1.10) 
This equation admits a static solution a that satisfies 

^-m^^a-Aa^^O. (1.11) 

Let us take the quartic coupling A to be small. In particular we assume A <^ w?^/ Jq. Then 
the static solution for a, to first order in A is 



^ - 7 i I ■ (1.12) 

H . 



The nonzero vacuum expectation values (VEV) for a and d^a imply that the scalars are in 
the condensate phase, with nonzero number density: Jq = doaa^ ^ 0. 

Prom the Lagrangian (1.8) we can calculate the spectrum and propagation of perturba- 
tions above the condensate. We find a heavy mode of mass 2mH which we ignore as it is 
beyond the scope of the low energy theory (since we assume that itih ^ •^o''^)) ^ well as a 
light mode. The dispersion relation for the light mode is as follows: 

a;2~^^ + ^^k2 (113) 

Here we have taken the limits A Jo <tC rriff and ^ m^. In the absence of the self- inter ation 
term, i.e. when A — * 0, this dispersion reduces to that for a free particle: u = k^/2mjy. 

The presence of the self-interaction term gives rise to the superfiuidity of the neutral gas. 
The long-wavelenth modes obey a linear dispersion relation: 




..,/^|k|. (1.14) 



The group velocity is then 



Such a linear dispersion relation (1.14) meets the Landau criterion for superfluidity [20]. At 
velocities less than Vg^ the system experiences no loss of energy due to its motion. Note that 
if the self-interaction coupling A were zero, then Vg^ would also be zero and there would be 
no velocity at which the system could sustain a non-dissipative flow. Self-interactions are 
essential to superfluidity. 

The internal energy, speciflc heat and other thermodynamic quantities of the interact- 
ing condensate also follow from the above dispersion relations. Using expression (1.14) as 
the energy e in the Bose- Einstein distribution (1.1), one can easily find the temperature 
dependence of the energy density of the gas of phonons {U oc T^) and of the specific heat 
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1.2.2 The non-relativist ic EFT 



The relativistic effective Lagrangian adopted in tlie previous section is not necessarily tfie 
most appropriate description of the low energy condensate of non-relativistic particles, al- 
though captures many of its significant features. It is overly restrictive in that it enforces 
Lorentz invariance- a symmetry we do not expect the low energy system to preserve. In this 
section we discuss a non-relativistic effective Lagrangian description of the neutral conden- 
sate. We will see that in this formalism the condensate retains the distinctive features found 
in the relativistic theory, namely the equivalent dispersion relation for the light mode. 

A non-relativistic effective order parameter Largangian that is consistent with the sym- 
metries of the physical system can be written as: 

Again, higher order terms can be included in this Lagrangian, however they are irrelevant 
for our discussions. The possible quadratic term, fx^R^*^, can be absorbed into the first 
term in (1.16) by redefinition of the phase of the scalar field. We switch notation, writing $ 
in terms of its modulus and phase $ = Sexp(ir). Written in terms of fields E and F, the 
effective Lagrangian (1.16) takes the following form: 

C^ff - -OoT S - ^— - ^ . (1.17) 

Varying w.r.t. E and F gives the following equations of motion: 

rriH 

\y;3 

V^E = 2m//EaoF + E(aF)2 + . (1.19) 

ruH 

Taking djV = 0, the first equation (1.18) gives the conservation of the scalar number density: 
9o(E^) = 0. Accordingly, we fix E^ = Jq. When A is nonzero, equation (1.19) has the 
following nonzero solution for 9oF: 

From the Lagrangian (1.17) we find the dispersion relation for the scalar perturbation 
above the condensate: 



This expression is identical to the approximate relation found for the light mode in the 
relativistic theory (1.13). 
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Chapter 2 

The Charged Condensate 



2.1 A description of charged condensation 

Consider a neutral system of a large number of nuclei each having charge Z, and neutralizing 
electrons. If the average inter-particle separation in this system is much smaller than the 
atomic scale, ~ 10^*^ cm, while being much larger than the nuclear scale, ~ 10^^^ cm, neither 
atomic nor nuclear effects will play a significant role. Moreover, the nuclei can also be treated 
as point-like particles. In what follows we focus on spin-0 nuclei with Z <S (hehum, carbon, 
oxygen), and consider the electron number-density in the interval Jo ~ (0.1 — 5 MeV)^. Then 
the electron Fermi energy will exceed the electron-electron and electron-nucleus Coulomb 
interaction energy. Moreover, at temperatures below ~ 10'^ K, which are of interest here, 
the system of electrons form a degenerate Fermi gas. 

Since the nuclei (we also call them ions below) are heavier than the electrons, the tem- 
perature at which they'll start to exhibit quantum properties will be lower. Let us define the 
"critical" temperature Tc, at which the de Broglie wavelengths of the ions begin to overlap 



where uih denotes the mass of the ion (the subscript "iJ" stands for heavy), and d denotes 
the average separation between the ions"*^. 

Somewhat below Tc quantum-mechanical uncertainties in the ion positions become greater 
than the average inter-ion separation. Hence the latter concept loses its meaning as a micro- 
scopic characteristic of the system; the ions enter a quantum-mechanical regime of indistin- 
guishability. Then, the many-body wavefunction of the spin-0 ions should be symmetrized. 
This unavoidably leads to a probabilistic "attraction" of the bosons to condense, i.e., to 

-'^The de Broglie wavelength above is defined as XdB = 27r/|k|, where k^/2mjj = 3fcBT/2. We define Tc as 
the temperature at which XdB — d. Note that this differs by a numerical factor of sJ^-k j'i from the standard 
definition of the thermal de Broglie wavelength, A = ^J1■K|mllT, that appears in the partition function of 
an ideal gas of number-density n in the dimensionless combination IS?n. 




(2.1) 
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occupy one and the same quantum state. We refer to the system of condensed nuclei and 

electrons charged condensate. 

In the condensate the scalars occupy a quantum state with zero momentum. Moreover, 
as we will show in section 2.3, small fluctuations of the bosonic sector have a mass gap, 
= (Ze^Jo/mj/)^/^, which exceeds Tc by more than an order of magnitude. Therefore, 
once the bosons are in the charged condensate, their phonons cannot be thermally excited. 
However, the gaplcss fcrmionic degrees of freedom near the fermi surface are thermally 
excited, and carry most of the entropy of the entire system [4-7,9]. 

For further discussion it is useful to rewrite the expression for Tc in terms of the mass 
density p = runJo measured in g/cm^: 

where the baryon number of an ion was assumed to equal twice the number of protons, 
A = 2Z (true for helium, carbon, oxygen). Thus, for p = 10® g/cm^ and helium-4 nuclei we 
get Tc ~ 10® K, while for the carbon nuclei with the same mass density we get Tc ~ 2 • 10^ 
K. 

The temperature at which the condensation phase transition takes place, T^ondens, is 
expect to be close to Tc but need not coincide with Tg. Calculation of Tcondens from the 
fundamental principles of this theory is difficult. However, we can obtain an interval in 
which Tcondens should cxist. For this we introduce the following parametrization: 

Tcondens — C -^c ) (^■^) 

where C is an unknown dimensionless parameter that should depend on density more mildly 
than Tc does. Numerically, this parameter should he in the interval 0.1 <^ C ^ 1^ the 
point C = 0.1 corresponds to the Bose- Einstein (BE) condensation temperature of a free gas 
for which T^^^^g^^ — LS/mnd^ is known from fundamental principles (see chapter 1). The 
condensation temperature in our system should be higher than T^^^^^^ since the repulsive 
interactions between bosons makes it easier for the condensation to take place [24]. In our 
case, these repulsive interactions are strong - the Coulomb energy is at least an order of 
magnitude greater than the thermal energy in the system. Hence, we would expect ( ^ 0.1. 
Moreover, we will show in section 2.3 that the quantum dynamics of the fcrmions screen the 
mass of the bosons, in effect lowering tuh- Thus the critical temperature defined by equation 
(2.1) can be greater due to this smaller effective mass. In what follows we will retain ( in 
our expressions, but use the value, C — 1) when it comes to numerical estimates. Perhaps 
a more accurate estimate of the condensation temperature can be obtained along the lines 
of [25]. 

The condensation will take place after gradual cooling only if Tcondens is greater than the 
temperature at which the substance would crystallize. A classical plasma crystallizes when 
the Coulomb energy becomes about ~ 180 times greater than the average thermal energy 
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(2.2) 



per particle [26-28]. This gives the following crystallization temperature^ 



T„y,t ^ (0.8 ■ 10=^Z^/3) K. (2.4) 

Note that the density dependence of is different from that of To-yst - for higher densities 
grows faster, making condensation more and more favorable! One can define the "equality" 
density for which Tcondens = T^yst: 

Pe,= (^)'^^°g/cm^ (2.5) 

For helium {Z = 2) we find peq — 10^ g/cm^, while for carbon {Z = 6) we find peq — 
10^ g/cm^ (as mentioned above, we use ( ~ 1). These results are very sensitive to the value 
of (] for instance, peq could be an order of magnitude higher if ( ^ 0.5. Regardless of this 
uncertainty, however, the obtained densities are in the ballpark of average densities present 
in helium-core white dwarf stars ~ 10^ g/cm^. For carbon dwarfs, they're closer to those 
expected in high density regions only [9]. 



T(°K) 




(a) Helium nuclei 



TTK) 




(b) Carbon nuclei 



Figure 2.1: Crystallization temperatures (blue, solid line) and condensation temperatures 
(red, dashed line) as a function of mass density for systems of helium and carbon nuclei. 

In Fig. 2.1 we plot the the crystallization temperature and condensation temperature 
as a function of mass density for systems of both helium and carbon nuclei, taking C = 1- 
The solid line indicates the crystallization temperature and the dashed line the condensation 
temperature. Typical core densities in white dwarf stars lie in the interval presented above, 
though the densities of most fall near the lower edge of this interval. We see that for 
helium-core white dwarfs the condensation temperature is significantly greater than the 
crystallization temperature for all given densities. The above considerations lead us to 
conclude that the dense system of the helium-4 nuclei and electrons may not solidify, but 
should condense instead. We will return to this system in chapter 3. 

^See chapter 3 for further discussion. 
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Fig. 2.1 also shows that for a carbon-core white dwarf, at typical densities, the critical 

temperature does not exceed the crystallization temperature and thus we would expect the 
nuclei to crystallize as usual. It can be seen that for superdense carbon white dwarfs with 
central density p ~ 10^° g/cur" the condensation temperature can exceed the crystallization 
temperature and thus these cores may also undergo condensation, however, this density 
is close to the neutronization threshold, as well as to the threshold where the relativistic 
gravitational instability would set in, so the existence of dwarfs with such a high average 
density is questionable. On the other hand, such densities may still exist in small regions in 
the very core of dwarf stars; some effects of this were studied in [9]. For oxygen-core white 
dwarfs (not shown) the crystallization temperature is always greater than the condensation 
temperature for relevant densities and thus we do not expect condensation to occur. 

Is the charged condensate the ground state of the system at hand? For the higher 
values of the density interval considered, the crystal would not exist due to strong zero-point 
oscillations. At lower densities, the crystalline state has lower free energy (at least near zero 
temperature) due to more favorable Coulomb binding. Hence, the condensate can only be a 
metastable state. The question arises whether after condensation at ~ Tcondens the system 
could transition at lower temperatures ~ Tcryst to the crystal state, as soon as the latter 
becomes available. 

In the condensate, the boson positions arc entirely uncertain while their momenta are 
equal to zero. In order for such a system to crystallize later on, each of the bosons should 
acquire the energy of the zero-point oscillations of the crystal ions. As long as this energy, ~ 
(Ze^ Jo/m//)^/^, is much greater than Tcryst, no thermal fluctuations can excite the condensed 
bosons to transition to the crystalline state. The latter condition is well-satisfied for all the 
densities considered in this work. There could, however, exist a spontaneous transition of a 
region of size Rc to the crystallized state via tunneling. The value of Rc, and the rate of this 
transition, will be determined, among other things, by tension of the interface between the 
condensate and crystal state, which is difficult to evaluate. However, for an estimate, the 
following qualitative arguments should suffice: the height of the barrier for each particle is 
(Ze^ Jo/mif)^/^ = m^, while the number of bosons in the Rc region ~ R^Jq/Z. Hence, the 
transition rate should scale as exp^—m^JoR'^/Z). Since we expect that Rc > the rate 

is strongly suppressed for the parameters at hand. 

2.2 The relativistic EFT 

To described the charged condensate in terms of an effective field theory, we start by consid- 
ering the simplest model that exhibits the main phenomenon: a generic, highly dense system 
of charged, massive scalars and oppositely charged fermions at zero temperature and infinite 
volume. The scalar field could be a fundamental field or a composite state in the regime 
that compositeness does not matter. The gauge field could be a photon or any other U{1) 
field. 

The scalar condensate is described by the order parameter 0. A nonzero vacuum expec- 
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tation value (VEV) of implies that the scalars are in the condensate phase. We adopt 
a relativistic Lorentz-invariant Lagrangian which contains the charged scalar field (f) with 
right-sign mass term mjj > 0, the gauge field yl^, and fermions "i/^ ^ with mass mp: 

>C = + I^M'/'r - mU*<l> + ^{i^D^ - mp)^ + fiF^''^ . (2.6) 

The covariant derivatives in (2.6) are defined as (9^ — ig(f,A^ for the scalars, and — ig^A^ 
for the fermions. Their respective charges, and gqr, are in general different. For simplicity 
we take g = g^ ^ -g^^. 

The chemical potential /x^ is introduced for the global fermion number carried by 
for example lepton or baryon number. The fermions in (2.6) obey the conventional Dirac 
equation with a nonzero chemical potential. In particular, a self- consistent solution of the 
equations of motion implies that 



I^F = ykj, + ml + gAo , (2.7) 

where kp denotes the Fermi momentum of the background fermion sea. The Fermi momen- 
tum is related to the number density of fermions Jo as follows: kp = (Stt^ Jq)^/^. A nonzero 
chemical potential implies a net fermion number density in the system. Since the fermions 
are electrically charged, they set a background electric charge density. Such charged fermions 
would repel each other. In our case, however, the fermionic charge will be compensated by 
the oppositely charged scalar condensate, as we will show below. At distance scales that are 
greater than the average separation between the fermions their spatial distribution can be 
assumed to be uniform. Then, the background charge density due to the fermions can be 
approximated as J/^ = 'ip^f^ip — Jo^f^o, where Jo is a constant, fixed by (2.7). One should take 
into consideration effects due to fermion fiuctuations above Jo. For now, however, we will 
assume that the fermions are frozen in "by hand" and thus the averaging procedure is valid. 
In later sections we will consider effects due to the dynamics of the fermions including their 
quantum loops. 

Because the system also has a conserved scalar current, we can associate with it a chemical 
potential Hs- For the Hamiltonian density, the inclusion of a chemical potential for the scalars 
results in the shift H ^ H' = H - /isJo'"^^^', where Jg^^'^''^'" = -i[{Do(f))* (f) - 0*(£'o0)] is the 
time component of the conserved scalar current. For the Lagrangian density this shift can 
be written as a shift in the covariant derivative for the scalar — > D'^ = — i/isSnQ. In 
what follows primed variables Tl', CJ will refer to those variables which include a nonzero 
chemical potential for the scalars. 

We have not included a quartic interaction term for the scalar A (0*0)^ in the Lagrangian 
(2.6). This term could be present, but it is straightforward to check that our results will not 
be affected as long as A Jo ^ m^. The VEV of the scalar is fixed by it being energetically 
favorable for the bulk of the condensate to be neutral. 

In general, the scalar field could have an additional Yukawa term, q{(^*'\\)]T't\)2 + h.c). 
where g is a coupling, F denotes either the 1 or 275 matrix depending on the spatial parity of 
0, and V'i_2 denote fermions with different ^7(1) charges that render the Yukawa term gauge 
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invariant. One or both of these fermions could be setting the background charge density Jq. 
A fermion condensate Jg'^ = {ipiip2 + h.c), if non-zero, could act as a source for the scalar. 
In order for this not to significantly change our results, the condition qJ^''^ ^ m\j{(f)) should 
be met.^ Due to this Yukawa coupling the scalar can decay. In order for the condensate 
phase to form in the first place, the "condensation time" must be shorter then the 

lifetime of the 0. 

In the case of a system of charged nuclei and electrons the Yukawa terms are forbidden 
by the gauge U{1) symmetry, so we will not be discussing them here. 

In addition, the effective Lagrangian (2.6) could contain the dimension-5 operator oc 
0*0^^. Such a term would renormalize the mass of the fermions. However, as long as 
the coupling constant that multiplies this term is sufficiently small this effect will not be 
significant. This is the case when the VEV of is much smaller than the UV cutoff scale by 
which the above dimension-5 operator is suppressed. The latter conditions is fulfilled in our 
case since Jo <^ m\. 

The complex order parameter can be written in terms of a modulus and a phase 
4> = -^(J e*". As per the discussion above, we treat the fermions as a fixed background density 
which couples to the gauge field as —gA^J'^. In terms of these variables the Lagrangian 
density becomes 

^' = -i^i + K^M^)' + IMm + f^sS,o - d,af a' - \m\ - gA.J'^ . (2.8) 

From this form of the Lagrangian it is evident that a nonzero expectation value for Aq, or 
doa, or a nonzero chemical potential /is can give rise to a tachyonic mass for the scalars. 
This effective tachyonic mass makes it possible for the scalar field to condense, as we shall 
now show. 

Varying the Lagrangian with respect to and a gives the following equations of motion: 

-d^'F^, = g{gA, + fiJ,o-d,a)a^ - gj,, (2.9) 
Oa = [{gA^ + fisSf,o-d^af -ml]a. (2.10) 

The Bianchi identity for the first equation in (2.9), can also be obtained by varying the 
action w.r.t. a. This gives the conservation of the scalar current: 

jscalar ^ [^^^^ ^ ^^^^ _ q^^^ ^ q_ (2.11) 

We can express the potential in terms of the gauge invariant variable = ^4^ -|- ^d^^a. 
For a constant charge density, J^^ = Jo^no, the theory admits a static solution: 




(2.12) 



When a acquires a VEV, the gauge symmetry is spontaneously broken. The mechanism 
of symmetry breaking for the charged condensate differs from the abelian Higgs model in 

^The Yukawa coupling would also lead to the new terms in the fermion mass matrix. Depending on the 
specific situation, this may or may not impose additional constraints. 
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that here the scalar field has a conventional positive-sign mass term. Instead of giving a 

tachyonic mass to the scalars by hand, a nonzero expectation value for Bq or a nonzero 
chemical potential fig act as a tachyonic mass term. In particular, when {qBq) + fJ^s = itt-h, 
the scalar field condenses. In the bulk of the condensate the scalar charge density exactly 
cancels the fermion charge density: J^^^'^'^^ = {{qBq) + Hs){cr'^) — Jq. 

2.3 Spectrum of perturbations 

The uniform fermion background sets a preferred Lorentz frame. We study the spectrum 
and propagation of perturbations in this background frame. For this we introduce small 
perturbations of gauge and scalar fields, and r, above their condensate values as follows: 

B^{x) = -{mH - l^s)5,j,o + b^{x) , a{x) = J — + t{x) . (2.13) 
g V 

In the quadratic approximation, the Lagrangian density for the perturbations reads 

A = + i( V)' + + '^rn^ruH h^r + ... (2.14) 

Here f^^, denotes the field strength for h^, and we have defined the following mass 

m, = g\f^. (2.15) 
V ruH 

We have dropped all the fermionic terms as well as the cubic and quartic interaction terms 
of fe's and r. This procedure is valid assuming that the perturbations are small compared 
to their condensate values, i.e. gbo <^ tuh and r -C \fj^fniH. The last term in (2.14) is 
Lorentz violating and is a consequence of having introduced the background fermion charge 
density. 

It is also worth noting that the chemical potential for the scalars /is has disappeared 
from the Lagrangian (2.14). In other words, the Lagrangian of perturbations is insensitive 
to whether the original theory (2.8) contained a nonzero chemical potential for the scalars, 
or whether a nonzero VEV of the gauge invariant field B^ is responsible for the condensation 
of the scalars. Both cases give rise to the same spectrum of small perturbations. 

Calculation of the spectrum of the theory is non-trivial but straightforward. We briefiy 
summarize the results. First, 6o is not a dynamical field, as it has no time derivatives in 
(2.14). Therefore, it can be integrated out through its equation of motion, leaving us with 
the equations for three polarizations of a massive vector hj (j = 1,2,3), and one scalar r. 
These constitute the four physical degrees of freedom of the theory. The transverse part of 
the vector hj obeys the free equation 

(□ + m2)6j = 0, where 6j = 6,- - ^(46^) . (2.16) 
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Therefore, the two states of the gauge field given by bj have mass rrij. Moreover, the 
frequency uj and the three-momentum vector k of these two states obey the conventional 
dispersion relation, cu^ = + k^. 

The longitudinal mode of the gauge field bj, and the scalar r, on the other hand, give 
rise to the following Lorentz- violating dispersion relations (valid for nu/ ^ 0) 

u;| = k^ + 2ml + ± ^4k2m2, + (2m|, - ^m^)^ . (2.17) 

The r.h.s. of (2.17) is positive, thus the condensate background is stable w.r.t. small per- 
turbations. Both of these modes have masses which can be obtained by putting k = 0. For 
one mode this mass coincides with m^. We refer to this mode as the longitudinal mode or 
the phonon, though in reality it is a linear combination of the scalar r and the longitudinal 
gauge boson The other mode has a mass vrtg = 2m h corresponding to the creation of a 
particle-antiparticle pair of scalars. We refer to this as the scalar mode. Interestingly, the 
group velocities of the transverse and longitudinal modes of the massive vector boson are 
different. For mn "rrij, and for an arbitrary k, the fastest ones are the transverse modes, 
they're followed by the scalar, and the longitudinal mode is the slowest. All group velocities 
are subluminal. 

In the non-relativistic limit, for mH ^ m^, the dispersion relations (2.17) simplify to: 

ujI ~ Am]j + 2k2 , ujl^mi+ \ , , (2.18) 

The first of these corresponds to the creation of a particle-antiparticle pair. The momentum 
term in this dispersion relation corresponds to the energy required for the propagation of one 

of these particles. The second dispersion relation is more unusual. To better understand it we 
consider the decoupling limit. When electromagnetic interactions are turned off, g ^ and 
thus the mass m^ — > 0. Accordingly the second dispersion relation becomes a;_ = \s? / (2mj^). 
This is exactly the energy required to propagate a free boson. Thus the m-y terms in this 
dispersion relation can be thought of as a consequences of the electromagnetic interactions. 
The mass gap m^ in the bosonic spectrum means that the Landau criterion is automatically 
satisfied and thus the bosons exhibit superfiuidity. The dispersion relations for uj- in (2.18) 
also exhibits roton-like behavior (more on this in the next section). 

In our discussions so far we have treated the fermions as a fixed charge background 
Jjji = <^^/iO• We relax this assumption now and introduce dynamics for the fermions via the 
Thomas-Fermi (TF) approximation. We consider the corrections to the spectrum of small 
perturbations due to these dynamics. 

The fermion number density is governed by the constant chemical potential ^p: 

liF = ^(37r2Jo(a;))2/3 + m| + gAo{x) . (2.19) 

Here we have related the local number density of fermions to the Fermi momentum via 
Jo(x) — kpixY /{Stt"^). In this way the number density of the fermions Jq gets related to the 
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electric potential Aq. For relativistic fermions 

Mx) = ^{fiF - gAo{x)f . (2.20) 

Consequently, fluctuations in Jq{x) can be expressed in terms of fluctuations in the potential 
bo{x). As a result, the coefficient in front of 6q gets modified as compared to (2.14) 

A = ~C + + >o&o - + 2m^^^ boT + ... (2.21) 

where ml = + {g'^/n'^) (Stt^Jq)^/^. The latter term is simply the Debye mass squared. 
The introduction of the fermion dynamics via the TF approximation breaks the degeneracy 
between the "electric" and "magnetic" masses of the gauge field. 

Calculation of the spectrum is once again straightforward. The two transverse compo- 
nents of the gauge field are not affected by the addition of the fermion dynamics. They still 
propagate according to the usual massive dispersion relation o;^ = + k^. The dispersion 
relations of the longitudinal and scalar modes (2.17) become 

± 4/4k2 ^+ — 2 ^ + ^ O 2 

y mg L mg 2 \ 2ml ) 

where M = ^m^m^. The r.h.s of (2.22) is positive for arbitrary k and thus the charged 
condensate background is stable w.r.t. small perturbations. All the group velocities obtained 
from (2.22) are subluminal. 

The solution with the minus subscript corresponds to the longitudinal component of the 
massive vector field, with a;^(k = 0) = m^. Though the dispersion relation is different with 
the introduction of the fermion dynamics, the mass gap for this mode remains the same. 
The solution with the plus subscript corresponds to the scalar mode and its mass squared 
in this frame is cc;^(k = 0) = 4M^/r7iQ. Prior to introducing the fermion fluctuations, the 
mass of this mode was given by oj\{k. = 0) = Arn^jj — 4M^/m^. As ml is greater than 
m^ by definition, we see that the effect of the fermion dynamics is to lower the effective 
mass of the bosonic particle-antiparticle pair. Hence the effective mass of the condensed 
bosons get screened due to fermion quantum effects, suggesting that the condensation phase 
transition temperature may actually be even higher than what we use here, and adopting 
C = 1 in chapter 1 may be a conservative choice. We cannot, however, directly use the 
above expression for the effective mass since the effective theory used for its derivation is 
not expected to be reliable at the scale of the mass itself. 

In the limit that the fermions are non-dynamical (i.e. are frozen "by hand" or some 
other dynamics), then mo m^, and the solutions reduce to the ones obtained in (2.17). 
However, for most physical setups we will find that the difference between mo and is 
greater than m^. Therefore, the fermion dynamics introduce an additional screening of the 
electrostatic interactions. In fact, this is just the usual Debye screening. 



(2.22) 
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2.4 Screening of electric charge 



As a next step we'd like to discuss the screening of electrically charged impurities in the 
condensate. To determine the screening length we consider a small, spherically symmet- 
ric object with a nonzero charge placed in the condensate. Outside of the charge, in the 
condensate, the equations of motion for a static Aq and a as derived from (2.9) and (2.10) 
are: 

-V'^ + g{gAo + ij,)a^ = ^ Jo , - VV = [{gAo + Hs? - mjj]a . (2.23) 

In terms of the perturbations of the fields above their condensate values (2.13), we can write 
these equations as 

V\ = m^fco + 2MV , - VV = 2M\ , (2.24) 

where again we have defined M = y/mHfn^. The above equations are valid as long as we 
focus on solutions that satisfy r <C \/ Jo /niH and gbo <^ mn- 

The regime of physical interest is the one in which uih ^ J^^^ ■ This will be apphcable 
to the system of helium-4 nuclei and electrons in helium white dwarf stars. In this case 
M 3> m^, and we can neglect the first term on the r.h.s. of the equation for in (2.24). For 
large r we require that 6o, t ^ 0. The boundary conditions select the decaying functions: 

hoir) ~ [cisin(Mr) + C2Cos(Mr)] , (2.25) 

r 

r(r) ~ [-ci cos(Mr) + C2 sin(Mr)] . (2.26) 

r 

The constants Ci and C2 are to be determined by matching these solutions to those in the 
interior of the small, charged object. Thus, for a probe particle, the screening occurs at scales 

1 /3 

greater than 1/M. When rriH ^ Jq , 1/M is shorter that the average inter-particle separa- 
tion d oc Jq^^^- Although this strong screening may well be a reason why the condensation 
of charged bosons takes place in the first place, this statement needs some qualifications. It 
may seem that the exponent exp(— Mr) is due to a state of mass M. The distance scale 1/M 
is shorter than the average inter-particle separation - an effective short-distance cutoff of the 
low-energy theory. Then, a state of mass M, if existed, would have been beyond the scope of 
the low-energy field theory description, and the above potential would have been unreliable. 
There is an explanation of this scale in terms of a cancellation between the potentials due 
to the two long-wavelength modes - Coulomb and "phonon" quasiparticles - both of which 
are much lighter than 1/d, and are well- within the validity of the effective field theory. So 
the obtained result is well within the scope of the effective field theory for r > d. 

To study this effect in more detail we start by calculating the gauge boson propagator. 
For now we once again treat the fermions as "frozen in:" = Jq^ixo- We will return to their 
dynamics at the end of this section. The gauge boson propagator can be determined from 
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the Lagrangian of small perturbations (2.14). It is useful to integrate out the r field. The 
remaining Lagrangian takes the form: 



(2.27) 



This Lagrangian contains four components of b^, and no other fields. The first two terms in 
(2.27) are those of a usual massive photon with three degrees of freedom. The last term is 
unusual, as it gives rise to the dynamics to the timelike component of the gauge field. This 
term emerged due to the mixing of bo with the dynamical field r in (2.14), and since we 
integrated out r, bo inherited its dynamics in a seemingly nonlocal way. 

This form of the Lagrangian (2.27) is useful for calculating the propagator. Indeed, the 
inverse of the quadratic operator that appears in (2.27) has poles which describe all the four 
propagating degrees of freedom. The full momentum-space propagator is given by 
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(2.28) 
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where p is the four-momentum and uj = po- In the limit that M 
propagator describes a usual massive vector boson: 
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For nonzero M the propagator is modified by Lorentz- violating terms. 

Sandwiched between two conserved currents and J^, the propagator takes the form: 
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(2.30) 



The poles of this propagator describe two transverse photons with mass m-y, one heavy mode 
with mass 2mij, and a light phonon with mass m-y. Their dispersion relations are those found 
in (2.16) and (2.17). 

In particular, we are interested in a static potential for a point source. This can be 
obtained from the propagator (2.30): 



y(k) =L'oo(k,a; = 0) 



1 



1 



k2 + 



k2 + m2 + k2(k2 + m2)2/4M4 ■ 



(2.31) 



The first term on the r.h.s. can be thought of as an repulsive screened Yukawa potential 
while the second term can be interpreted as an attractive potential due to a phonon. This 
second term in (2.31) has three poles. The residue of one pole exactly cancels that of the 
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first term of (2.31). The remaining two poles describe both the heavy state of mass 2mH 
which is unimportant for the low-cncrgy physics, and a hght state, which actually is the 
phonon. It's the light mode found in (2.17) that belongs to the spectrum of the low-energy 
effective field theory. For simplicity of the discussions, we'll be using a somewhat imprecise 
language by calling the whole second term in (2.31) the phonon contribution. 

The phonon in this case is a collective excitation of the motion of charged scalars within 
the fermion background. The cancellation due to this light mode gives rise to the exponential 
exp(— Mr), and not a hypothetical state of mass M. At scales larger than 1/M, which are of 
primary interest, the phonon potential cancels the gauge potential with a high accuracy. This 
cancellation is reliable at scales that are much greater than Jq^^^ ^ M~^, and takes place 
already at scales that are much shorter that the photon Compton wavelength ^ Jg 

In a Lorentz-invariant theory having a negative sign in front of a propagator, such as the 
one in the second term of (2.31), would suggest the presence of a ghost-like state. However, 
this is not the case in a Lorentz- violating theory described by our Lagrangian (2.14). The 
fact that there are no pathologies in (2.14), such as ghost and/or tachyons, can be seen by 
calculating the Hamiltonian density: 

^2 = + y + \p'r + l^d^rf + \mt,h] + 2^(a,7r, - 2m^m,r)^ . (2.32) 

Here, tTj = —foj and = OqT. The Hamiltonian is positive semi-definite. Hence, no ghosts 
or tachyons are present. Moreover, consistent with ones expectation, the second term in 
(2.31) disappears in the hmit M — > 0, where the Lorentz invariance of (2.14) is restored. 
The static potential (2.31) can be simphfied: 

V^(k)=(k^ + m^+^)'\ (2.33) 

The first, second, and third terms on the r.h.s. of (2.33) are due to the respective terms in 
(2.27). Interestingly, when mn ^ fn^, there is no scale at which the photon mass term in 
(2.33) would dominate: for > the mass term is sub-dominant to the k^ term, while for 
k^ < it is sub-dominant to the Af^/k^ term. The term M^/k^ is coming from the phonon 
cancellation and gives rise to significant modification of the propagator in the infrared. 

The above described properties of the propagator can be seen by calculating the coordi- 
nate space potential from (2.33). We find, in agreement with (2.25): 



V{r) ^ {Q,eQ,e) j ^6*^^(7(0; = 0, k) oc ^i^H^^^^cos(Mr) , (2.34) 

in which we assumed that r = |x| ^ 1/Af and ^ Af. This potential is sign-indefinite 
and undergoes modulated oscillations between repulsion and attraction. There are an infinite 
number of points in the position space where the force between classical charges would vanish. 
These are points where 

dV 

— (r = r„) = 0, n = l,2,... (2.35) 
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Any two charged probe particles separated by a distance Jq , where our calculations 
are reliable, would stay in a static equilibrium as long as V^(r„) < 0. However, the potential 
is too shallow and any realistic temperature effects would kick the probe charges out the 
potential wells in (2.35). 

Furthermore, the finite temperature corrections could modify the properties of the con- 
densate itself, however, in this particular case, due to a high mass gap, the main prop- 
erties of the condensate should remain valid at temperatures well-below the condensa- 
tion point. For instance, in white dwarfs with temperature 10^ — 10'' K we would ex- 
pect the dominant temperature-dependent corrections to the potential to be proportional to 
T/Jq/^ ~ (10-^ - 10-^) < 1, which are negligible. 

Before proceeding to the next section we make two comments. First, in the uih oo 
limit one would expect the heavy scalars to decouple. It is not exactly clear from (2.31) how 
such a decoupling takes place, and what is its interpretation. In the limit rriH — oo, which 
implies M — > oo, the phonon effects should disappear. This certainly is the case in the full 
amplitude discussed before. However, taking this limit in (2.31) (or in (2.33)) results in a 
vanishing of the whole potential. This is an artifact of using the static approximation and 
can be understood in the following way: the phonon mixes with the timelike component of 
the gauge field, and because of this acquires an instantaneous part. Then, the instantaneous 
parts in (2.31) cancel between the gauge and photon contributions. However, the dynamical 
part of the phonon also reduces to zero, as the group velocity of the phonon vanishes in the 
niH — oo limit. 

This can be seen by looking at the dispersion relation for the phonon which was obtained 
in (2.18). For the relevant momenta ^ the dispersion relation for a;_ gives the phonon 
group velocity 



This vanishes in the rriH — > oo limit. 

Note that for ~ itl'^/2 the phonon group velocity also vanishes for finite itlh- This 
describes a state of a nonzero momentum but zero group velocity. The energy of this state 
is also nonzero, and to a good approximation equals to m^. These properties are similar to 
those of a roton in superfiuid hehum II. Moreover, for excitations with > m^/2 the group 
velocity is positive, while in the opposite case, < m^/2, it becomes negative (i.e. the 
direction of the momentum and that of group velocity are opposite to each other). These 
excitations resemble the positive and negative group velocity rotons in superfiuid helium II. 

The second comment concerns the limits of applicability of the linearized approximation. 
The expansion in (2.27) is only valid when the perturbations r and bo are much smaller 
than their condensate values. For scalar field this means that r <g yfj^frn^ and for the 
gauge field gb^^ ^ ^u- In the limit van ~^ oo, the domain of applicability of the linearized 
results shrinks to zero. This suggests that the geometric size of the region in which one can 
meaningfully talk about the charged condensate should be greater than a certain critical size 
that scales as ~ {\/~JoJmH)~^ ■ The latter tends to infinity as — > oo. 




(2.36) 
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In our treatment of the screening of electric charge we have treated the fermions as 
"frozen in." Such an approximation would be physically justifiable if, say, the fermions were 
fixed in a crystal lattice. However, in many physical circumstances, including the helium 
white dwarf system to be discussed later, this is not a good approximation. The fermion 
fluctuations should be taken into account. This was done in the previous section for the 
spectrum of small perturbations using the Thomas-Fermi (TF) approximation. The result 
was that the timelike component of the gauge field acquired an contribution to its mass. 
This mass term modifies the first of the equations of motion for small perturbations (2.24): 

V^6o = mlbo + 2MV , - VV = 2M\ , (2.37) 

where ml is defined as above as the sum of the photon mass squared and the Debye mass 
squared. The regime of physical interest, uih ^ Jq^^-, corresponds to M ^ mo. If we 
again calculate the potential outside a static probe charge using these modified equations, 
the rriQ term in the above expression for 60 is subdominant compared to the M term. Thus 
the potential obtained in (2.25) for a probe charge is still valid with the inclusion of this 
additional mass term. 

However, the TF approximation does not capture the significant property of the fermion 
system related to the possibility of exciting gapless modes near the Fermi surface. We 
can incorporate these effects into our results by calculating the one-loop correction to the 
propagator (2.28). For this, we restore back in the Lagrangian (2.27) the fermion kinetic, 
mass and chemical potential terms and, upon calculating the gauge boson propagator, we 
will take into account the known one-loop gauge boson polarization diagram. This diagram 
is suppressed by an additional power of the electromagnetic coupling constant ctem = e'^/in, 
and one would expect the quantum correction to be insignificant. However, this is not the 
case for the following subtle reason. The one-loop correction introduces branch cuts in the 
propagator, which give rise to additional contributions to the static potential in the position 
space. These additional terms have oscillatory nature with a power-like decaying envelope. 
Even though they are formally suppressed by 0{al), to a good approximation they end up 
being 0{1), and can dominate over the exponentially suppressed term at sufficiently large 
distances. 

The static potential obtained from the {00} component of the propagator Z^oo with the 
one-loop correction is given by 

V{k) = Doo{u;^0,k)^ lk^ + m''^ + -^ + F{e,kF,mF)] . (2.38) 

The function F{k'^^ kp, mp) is due to the one-loop photon polarization diagram and includes 
both the vacuum and fermion matter contributions {kp again denotes the Fermi momentum). 
A complete expression for F{k'^,kp,mp) can be found in Ref. [29]. We concentrate on the 
expression for F{k^,kF,mp) in the massless {mp — 0) limit that is a good approximation 
for ultra-relativistic fermions: 

164 + Ml^zM,„(a^)^_Pl^^ll_p:).V (2.39) 
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Here /Xq stands for the normalization point that appears in the one-loop vacuum polarization 
diagram calculation. The function F introduces a shift of the pole in the propagator, corre- 
sponding to the "electric mass" of the photon. This part of the pole can be incorporated via 
the TF approximation, as was done above. In addition, however, the function F also gives 
rise to branch cuts in the complex |k| plane (see [30] for the list of earlier references on this). 

In analogy with the propagator found above (2.31), we can decompose the static potential 
as follows: 

^(k,a; = 0) ^ (k^ + + F)-^ - (k^ + + F + ' ■ (2.40) 

The first term in (2.40) is just the instantaneous screened-Coulomb (Yukawa) potential of a 
massive photon with the one-loop polarization correction. Our main interest is at distances 
smaller than A sphere of radius encloses many particles within its volume since 

» Jq^^^. At these scales, the first term in (2.40) can be approximated by: 



7 

1 



k2 + F 



(2.41) 



The above expression has a regular pole corresponding to the acquired "electric" mass of the 
photon due to the polarization diagram. The contribution of this pole would give rise to an 
exponentially decaying potential e~"*«i''/r, where iriei ^ efj^j. This is just an ordinary Debye 
screening. 

However, as was mentioned above, the expression (2.41) also has branch cuts in the 
complex |k| plane for k = ±2kF. These branch cuts give rise to the additional terms in 
the static potential which are not exponentially suppressed, but instead have an oscillatory 
behavior with a power-like decaying envelope. In a non-relativistic theory they're known as 
the Priedel oscillations [30]. In the relativistic theory they were calculated in Refs. [29,31] 
(we follow here [29] and for simplicity ignore the running of the coupling constant due to the 
vacuum loop): 

^y^QiQ2<.sM2k^ (2 42) 

These branch cuts have a physical interpretation: since there is no mass gap in the fermion 
spectrum, a photon can produce a near-the-Fermi-surface particle-hole pair of an arbitrarily 
small energy and a momentum close to ±2/^^. The imaginary part of the one-loop photon 
polarization diagram should include the continuum of such near-the-Fermi-surface pairs. 
These are reflected as logarithmic branch cuts in the expression for F. 

Thus, if the phonon term (the second term) on the r.h.s. of (2.40) were absent one would 
have a power-like behavior (2.42) of the static potential at scales r < m~^. The phonon term, 
however, significantly reduces the strength of this potential. The result for this term can 
be calculated by directly taking the Fourier transform of (2.38). The dominant contribution 
comes from the branch cuts at A; = ±2kp. Drawing the contours around these cuts in 
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the upper half plane of complex |k| [29,30], one deduces the result. In the approximation 
M ^ kp ':$> m-y, which is relevant for our system, a static potential between like charges 
scales as 

4Qi<g2Qiem 4 sm{2kFr) 

The potential (2.43) is a generalization of the Friedel potential to the case when in addition 
to the fermionic excitations there are also collective modes due to the charged condensate. As 
this contribution to the overall potential is a result of a subtraction between the conventional 
Friedel term and the long-range oscillating term due to a phonon, its magnitude is suppressed 
by a factor of 16(A;ir/M)^, as compared to what it would have been in a theory without the 
condensed charged bosons (see [30] for the discussion of the conventional Friedel potential, 
and Ref. [19] for its recent detailed study in the presence of the charged condensate at finite 
temperature.)'* 

Nevertheless, AV in (2.43) dominates over the exponentially suppressed part of the to- 
tal potential found in (2.34), for separations between probe particles large enough for the 
effective field theory description to be applicable. The net static interaction in the charged 
condensate, set by (2.43), is very weak. It is, however, still much stronger than gravitational 
interaction between a pair of light nuclei. Although formally AV in (2.43) is proportional 
to ttg^, to a good approximation it is independent of ctem since oc g^imnJoY- 

The net potential takes the form 

Ae"^'" X 4acm A;isin(2A;i.r) \ 

Vstat = «e.Qig2 [-^cos{Mr) + J" ' • (2.44) 

The first, exponentially suppressed term modulated by a periodic function, is due to the 
cancellation between the screened Coulomb potential and that of a phonon [7]. The existence 
of such a potential due to cancellation between the photon and phonon exchanges was first 
pointed out in Ref. [16], in the context of superconductivity. 

More important, however, is the second term in (2.44) that has a long-range [7]. It 
dominates over the exponentially suppressed term in (2.44) for scales of physical interest, 
and exhibits the power-like behavior modulated by a periodic function. 

This potential (2.44) is not sign-definite. In particular, it can give rise to attraction 
between like charges; this attraction is due to collective excitations of both fermionic and 
bosonic degrees of freedom. This represents a generalization of the Kohn-Luttinger effect [33] 
to the case where, on top of the fermionic excitations, the collective modes of the charged 
condensate also contribute. 

In the charged condensate Cooper pairs of electrons can also be formed. However, the 
corresponding transition temperature and the magnitude of the gap, are suppressed by a 
factor of exp(— l/Cg^j), where e^^f is proportional to the value of the inter-electron potential 
that contains both screened Coulomb and phonon exchange. The fact that this potential 

^Note that for spin-dependent interactions the same effects of the charged condensate would give a 
generalization of the Ruderman-Kittel-Kasuya-Yosida (RKKY) potential [32]. 
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(2.43) 



has an attractive domain, though very small, can be seen from the static potential found 
above (2.44); the latter is suppressed by a power of a large scale M. Furthermore, taking 
into account the frequency dependence of the propagator in the Eliashberg equation does not 
seem to change qualitatively the conclusion of a strong suppression of the Green's function 
and pairing temperature. 

Hence, even though the bosonic sector (condensed nuclei) is superconducting at rea- 
sonably high temperatures < 10^ interactions with gapless fermions could dissipate the 
superconducting currents. Only at extremely low temperatures, exponentially close to the 
absolute zero, could the electrons also form a gap leading to the superconductivity of the 
whole system. In the present work we consider temperatures at which electrons are not 
condensed into Cooper pairs, and ignore the finite temperature effects. 

Finally we note that the magnetic interactions are not screened at the scale Af Instead, 
as is clear from the Lagrangian (2.21), the magnetic interactions are screened at scale of the 
Compton wavelength of the massive photon. This scale will end up being much greater that 
the average inter-particle separation d. 

This may seem somewhat puzzling since the one-loop fermion correction to the transverse 
part of the photon propagator may be expected to introduce corrections that are of the order 
(ttIq — m^), which would dominate over any effect of the order m^. This would be the case, 
for instance, in a plasma where the fermion loop would determine the plasma frequency. 
However, in the charged condensate the issue is more subtle, as was already seen in the 
zero-zero component of the photon propagator: the modification of the propagator due 
to the condensate does not simply reduce to a shift of the pole by m^, but instead gives 
rise to an additional infrared- sensitive momentum- dependent term in the propagator. This 
modification is such that the fermion loop correction to the real part of the pole is negligible 
in comparison with the contribution due to the charged condensate. In other words, the 
fermion-loop correction to the mass of the the longitudinal mode is negligible (the correction 
is of the order of ml/rri^ -C 1). The transverse mode has to have a mass equal to that of 
the longitudinal mode, which is determined by m^, since there is no difference between the 
transverse and longitudinal modes at zero momentum. 
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Chapter 3 

Helium White Dwarf Stars 



3.1 Condensation versus crystallization 

White dwarf stars represent a final evolutionary state of low mass stars. Because white dwarfs 
have exhausted their thermonuclear fuel, they are no longer supported against collapse by 
the heat generated by fusion. Instead they are stabilized by the degeneracy pressure of the 
electrons balancing against the gravitational attraction of the ions. As a result white dwarf 
stars are very dense - they are roughly of the size of the Earth and their mass is on the order 
of a solar mass. Their central mass densities range from ~ (10^ — 10^) g/cm^, with most 
of them falling near the lower edge of this interval. Their cores consist of a neutral system 
of electrons and nuclei (ions), the interparticle separation between nuclei being in between 
the atomic scale (Angstrom ~ 10~^ cm) and the nuclear scale (Fermi ~ 10~^^ cm). Thus 
the electrons and the nuclei are unable to form neutral atoms, yet nuclear effects can be 
considered insignificant. 

What, then, is the state of matter in the cores of these stars? To a certain extent, the 
answer is known - it depends on the temperature T. At high temperatures the equilibrium 
state is a plasma of negatively charged electrons and positively charged nuclei. However, at 
lower temperatures the system undergoes significant changes. Generally, as the star cools 
below a critical temperature Tcryst, the plasma becomes strongly coupled enough for the 
ions to crystallize [2]. This is the case for the majority of white dwarfs whose cores are 
composed of carbon or oxygen nuclei. However, for certain systems, in particular those 
whose cores are composed of helium nuclei, the de Broglie wavelengths of the ions begin to 
overlap significantly before the crystallization temperature is reached. In this case quantum 
effects can prevent the crystallization transition. Instead, the interior of the helium dwarfs 
may form a quantum state in which the electrons form a degenerate Fermi liquid and the 
charged helium-4 nuclei condense into a macroscopic state of large occupation number - the 
charged condensate. 

To show this, let us first consider the cooling process in the core of a white dwarf while 
ignoring quantum mechanical effects for the nuclei. We will treat the electrons as a degenerate 
Fermi liquid, as the Coulomb energy of a pair of electrons in a white dwarf is smaller than 
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the Fermi energy. For a typical dwarf star, below a certain temperature Tcrysu the Coulomb 
energy of a pair of ions will significantly exceed their classical thermal energy. In order to 
minimize their Coulomb energy, the ions will arrange themselves into a crystal lattice [2]. 
The crystallization temperature is characterized by the dimensionless ratio of the average 
Coulomb energy of a pair of ions to the thermal energy (see, e.g., [34]) 

p _ Ecoulomb _ {ZeY 1^ 



3 



Erhermal 47r(i T 



Here e denotes the electric charge, Ze is the charge of a nucleus, and we have set the 
Boltzmann constant ks = 1- The interparticle separation of the nuclei is given hy d = 
(47rJo/3Z)~^/^ where Jo is the electron number density. The second equality in (3.1) assumes 
the validity of the classical approximation. 

Numerical studies have indicated that when the temperature drops low enough so that 
r > 180, the ion plasma becomes strongly coupled enough for the system to crystallize (for 
earlier works see Refs. [2,26], for later studies see [27,28,35] and references therein). This has 
direct relevance to white dwarf stars: it is expected that in most white dwarfs, consisting of 
carbon, oxygen, or heavier elements, the crystallization transition takes place in the process 
of cooling (for a review, see, e.g., Ref. [3]). 

The above discussions were classical. As the star cools, quantum effects may become sig- 
nificant before crystallization sets in. For instance, in certain systems the zero-point energy 
of the ions can exceed the classical thermal energy T before the crystallization temperature 
is reached. The Debye temperature 0d gives the scale at which these considerations become 
relevant: 

/ T \ 1/2 

e^ = Qp, Qp = Ze(^-^j , (3.2) 

where Q,p is the ion plasma frequency and mn is the ion mass. The plasma frequency is related 
to the zero-point energy of the ions by ojq = Qp/y/S. Often, Oo may significantly exceed the 
crystallization temperature Tcryst [36]. In such cases, quantum zero-point oscillations should 
be taken into account. This seems to delay the formation of quantum crystal, lowering Tcryst 
from its classical value at most by about ~ 10% [35]. Since this is a small change, we will 
ignore it in our estimates. 

There is another scale at which quantum effects must be taken in consideration. For 
some systems, in particular those with light ions or those that are extremely dense, the de 
Broglie wavelengths of the ions will begin to overlap before the crystallization temperature 
is reached. The de Broghe wavelength is given by AdB = 27r/|k| where \s.^/2mH — ST/2. The 
overlap of the de Broglie wavelengths indicates that the quantum mechanical uncertainty in 
the position of the particles is greater than the average inter-particle separation. The critical 
temperature Tc at which this overlap occurs is given by 

47r2 ^ ^ 
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Below Tc the ions must be treated as indistinguishable particles obeying Bose-Einstein statis- 
tics. The statistics of indistinguishable bosons at low temperatures lead to the condensation 
of the bosons into the lowest energy quantum state. We discussed this in detail in chapter 1. 
Thus for systems in which Tc > Tcryst, instead of crystallizing, the system may condense into 
a macroscopic zero-momentum quantum state with large occupation number - the charged 
condensate.^ 

Dynamically, the condensation proceeds in the conventional manner: at temperatures 
higher than the condensation temperature most of the states are in thermal modes, and 
HsiT) is less than mn- As the temperature drops, HsiT) increases. (Here we ignore the 
temperature dependence of mu which would be present due to the quantum loop effects). 
As T — > Tc a- significant fraction of the modes ends up in the zero momentum ground state, 
while HsiTc) asymptotes to ttih- 

We wish to consider a system in which the quantum effects described above are maximal. 
Because helium-4 nuclei are lighter and have lower charge than carbon and oxygen nuclei, 
we expect them to be more sensitive to quantum effects. While the majority of white dwarf 
stars have cores composed of carbon or oxygen nuclei, helium-core white dwarfs constitute 
a small sub-class of dwarf stars (see, [10,37] or references therein). Most of helium dwarfs 
are believed to be formed in binary systems, where the removal of the envelope off the dwarf 
progenitor red giant by its binary companion happened before helium ignition, producing 
a remnant that evolves to a white dwarf with a helium core. Helium dwarf masses range 
from ~ 0.5 M© down to as low as (0.18 — 0.19) M©, while their envelopes are mainly 
composed of hydrogen. The system is long-lived as the helium-4 nuclei are stable w.r.t. 
fission. Furthermore, some nuclear reactions that could contaminate the helium-4 cores by 
their products are suppressed. One of these is the neutronization process due to inverse 
beta-decay. In our case the electrons are not energetic enough to reach the neutronization 
threshold of the helium-4 nucleus, which is about 20 MeV. Moreover, we would expect that 
the rate of helium fusion via the triple alpha-particle reaction is suppressed as the so-called 
pycnonuclear reaction (i.e., nuclear fusion due to zero-point oscillations in a high density 
environment) rates are exponentially small [34]. Hence, the cores of these white dwarfs are 
expected to be dominated by helium-4 for very long time after their formation. 

Table 3.1 summarizes the relevant physical quantities for a helium white dwarf star with 
a typical mass density p ~ 10^ g/cm^. We see that the condensation temperature 10^ K is 
significantly greater than the crystallization temperature 3 ■ 10^ K. Thus we expect that the 
system of the helium-4 nuclei and electrons may not solidify, but may condense instead. For 
hehum-core white dwarfs with even denser cores, the effect is greater still. 

3.2 The non-relativistic EFT 

As we discussed in chapter 1, for a large number of quanta that form a macroscopic state, 
such as the heliuni nuclei in the charged condensate, the state can be described in terms of 

^Following arguments of the previous chapter, we take the condensation temperature to be roughly that 
of the critical temperature: Tcondens — Tc- 
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Table 3.1: Typical physical quantities for a helium white dwarf 



Physical quantity 


Numerical value 


Mass density 
Electron number density 
Separation between atomic nuclei 

Debye temperature 
Critical temperature 
Crystallization temperature 


10^ g/cm^ 
(0.13 MeV)=^ 
10^ fm 
4 ■ lO'' K 

10^ K 
3 • 10^ K 



an effective field theory of the order parameter, and its long wavelength fluctuations. The 
effective field theory should be constructed based on the fundamental symmetries of the 
physical system and the properties of the interactions involved. The relativistic effective 
Lagrangian adopted in the previous chapter (2.8) is not the most appropriate description 
for the charged condensate in the cores of helium white dwarf stars, although captures 
many of its significant features. It is overly restrictive in that it enforces Lorentz invariance 
- a symmetry we do not expect the low energy system of helium-4 nuclei and electrons 
to preserve. In addition it contains a heavy mode of mass which we would expect to be 
beyond the scope of a low energy theory. In this section we discuss an low energy, non- 
relativistic effective Lagrangian description of charged condensation, and in particular its 
application to the system of helium-4 nuclei and electrons described above. We will see 
that in this formalism the charged condensate retains the distinctive features found in the 
previous chapter, namely equivalent Lorentz- violating dispersion relations for the massive 
photon and the same strong screening of electric charge. 

We focus on the zero-temperature limit, even though realistic temperatures in helium 
white dwarfs are well above zero. The validity of the zero-temperature approximation is jus- 
tified a posteriori as follows: the spin-0 nuclei undergo condensation to the zero-momentum 
state; while they do so they cannot excite their own phonons since the latter are gapped 
with the magnitude of the gap being greater than the condensation temperature. On the 
other hand, the condensing charged bosons can and will excite thermal fiuctuations in the 
fermionic sector that is gapless. Therefore, all the thermal fiuctuations will end up being 
stored in the fermionic quasiparticles near the Fermi surface. For the latter, however, the 
finite temperature effects aren't significant since their Fermi energy is so much higher than 
the temperature, T j J^^ <S 10~^. We note that the finite temperature effects, in a general 
setup with condensed bosons, were calculated in Refs. [18, 19]. 

The electron/nuclei system in the cores of the white dwarf has three relevant mass scales: 
the mass of a nucleus mui the electron mass me and the electron chemical potential [ip. The 
mass of the nuclei is significantly greater than the other two mass scales. For the effective 
Lagrangian construction we consider scales that are well below the heavy mass scale m^f, 
but somewhat above the scale set by rn^, and ixp. Hence the electrons are described by the 
Dirac Lagrangian, while for the description of the nuclei we will use a charged scalar order 
parameter $(a;). 
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The effective Lagrangian for $ must satisfy the following requirements: (i) it must be 
consistent with the symmetries of the physical system: translational, rotational, and Galilean 
symmetries, a global U{1) symmetry for the conserved scalar number, and a local gauge in- 
variance; (ii) it should reproduce the standard Schrodinger equation for the order parameter 
in lowest order in the fields; (iii) it should obey an algebraic relation between the conserved 
current density and the momentum density: J, = {e/mH)Toj; (iv) it should give an appro- 
priate spectrum of Nambu-Goldstonc bosons in the decoupling limit. Such a Lagrangian was 
first proposed by Greiter, Wilczek and Witten (GWW) [38] in a context of superconductivity: 

A// = V Q(rDo$ - {Do^m - , (3.4) 

where Do = (do — i2eAo), Dj = (dj — i2eAj). The function P(x) stands for a general 
polynomial function of its argument. The coefficients of this polynomial are dimensionful 
numbers that are inversely proportional to powers of a short-distance cutoff of the effective 
field theory, V{x) = Y.n=o^nX''- 

The Yukawa terms are forbidden by symmetries in this case. However, in general one 
should also add to the Lagrangian terms unr^*^, A($*<l>)^/m|^, 

— 1/3 

Ai($*$)^^/(mj^ Jq ), and other higher dimensional operators that are consistent with all 
the symmetries and conditions that lead to (3.4) (the Yukawa term is not). Here unr denotes 
a non-relativistic chemical potential for the scalars. These terms are not important for the 
low-temperature spectrum of small perturbations we're interested in, as long as A, Ai < 1 
and Jo -C m^. However, near the phase transition point it is the sign of /inr that would 
distinguish between the broken and symmetric phases, so these terms should be included 
for the discussion of the symmetry restoration. We also note that the scalar part of (3.4) is 
somewhat similar to the Ginzburg-Landau (GL) Lagrangian for superconductivity. However, 
there are significant difi^erences between them, one such difference being that the coherence 
length in the GL theory is many orders of magnitude greater than the average interelectron 
separation, while in the present case, the "size of the scalar" $ is smaller that the average 
interparticle distance. 

The GWW effective Lagrangian can also describe charged condensation, given the appro- 
priate rescaling and reintcrprctation of the parameters of the theory. In the condensate where 
the VEV of $ is nonzero, we can express $ in term of a modulus and phase: $ = S exp(ir). 
Written in terms of fields E and F, the effective Lagrangian (3.4) takes the following form: 

= V (^{2eAo - dom' - ^(V,S)^ - ^(2eA, - a,r)^E^^ . (3.5) 

For the fermions we adopt the usual Dirac Lagrangian with a relativistic chemical potential 

Cf = ^{iYD^ - mF)i^ + liFi^H , (3.6) 

where the covariant derivative for the electrons is defined as — dfj, + ieAfj_. The elec- 
tron chemical potential //p is equal to the Fermi energy at zero temperature /if — — 
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[(Stt^Jq)^''^ + rripY^'^. There exists a static solution that follows from Lagrangians (3.5) and 
(3.6). In the unitary gauge F = 0, this solution has the form 

2eE2 = e Jo , = , V'(0)^1. (3.7) 

Since on the solution the argument of (3.5) is zero, the condition V'{0) = 1 is satisfied by any 
polynomial function for which the first coefficient is normalized to one: V{x) = x + C2X^ + .... 
The above solution describes a neutral system in which the helium-4 charge density 2eE^ 
exactly cancels the electron charge density —eJo- 

Let us now turn to the issue of fluctuations about the classical solution. We express E 
in terms of a perturbation r above the condensate value: 

E(x)=V^(|y^ + r(x)j . (3.8) 

The Lagrangian, including the gauge fleld kinetic term, expanded to second order in fields 
becomes: 

A = -Ji^,l-^(5,r)^ + ^C^2m^^o^ + ^(37rVo)^/^)A^ (3.9) 

—m'^A'^j + 2mHm^Aor . 

Here we have included fiuctuations in the electron number density via the TF approximation. 
This gives rise to the second mass-like term, i.e. the Debye mass, in front of Aq. The photon 
mass is defined as m-y = 2e Jo/2mH- 

We can compare this Lagrangian to the one obtained in the relativistic theory (2.21). 
(Here the perturbation 6^ is equal to A^, as for the gauge field we are expanding above the 
background A^ = 0.) The first difference between the two Lagrangians is the lack of a time 
derivative for r in the non-relativistic case. This is because the heavy field is decoupled 
in the low energy limit. This, obviously, does not change the conclusions about the static 
potential discussed in the previous chapter. 

The second difference between (3.9) and (2.21) is the presence of the expansion coefficient 
C2 in the former. This is due to the fact that the non-relativistic theory for the scalar sector 
is not required to be Lorentz invariant. While the mass term due to the electron fiuctuations 
also breaks the degeneracy between the "electric" and "magnetic" masses of the gauge field, 
this is because the electron number density sets a preferred Lorentz frame. In the absence of 
this term, a Lorentz invariant Lagrangian for the scalars would have to have C2 = l/{mHJci) 
and thus the "electric" and "magnetic" masses would be equal. Within the non-relativistic 
effective Lagrangian approach, however, the value of C2 cannot be fixed a priori. It must be 
fixed instead by the physics of the particular system being described. 

As long as C2 is chosen so that the "electric" mass of the photon does not exceed the scale 
M, then our results will not greatly depend on the particular choice of C2. To calculate the 
spectrum of perturbations for the non-relativistic theory, we set by hand C2 — l/{mHJo) for 
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simplicity. Then the Lagrangian for small perturbations above the condensate in the non- 
rclativistic theory is identical to the Lagrangian for small perturbations in the relativistic 
theory up to a time derivative for r: 

A = -\fI - \{d,Tf + \mlAl \m\A] + ^mum^A^r , (3.10) 

where again rriQ = m^+e^/n^ (Stt^ Jq)^/^. Calculating the spectrum of perturbations from this 
Lagrangian gives two transverse components of the gauge field which propagate according 
to the usual massive dispersion relation o;^ = + k^. In addition there is a longitudinal 
mode with the dispersion relation: 

For the last equality we have taken the low momentum limit. This dispersion relation should 
be compared to those found for the relativistic theory in the rriH ^ itl^ limit (2.18). We 
find that the heavy mode found in (2.18) is no longer present in the low energy theory, 
as we would expect. In addition, the term proportional to — m^k^/4m|f that appears in 
the longitudinal dispersion relation for c<j_ in (2.18) does not appear in the non-relativistic 
theory, since it is of a higher order in I/itlh- It was this term that was responsible for the 
roton-like behavior of the small fluctuations discussed in the previous section. 

Nevertheless, the longitudinal excitation exhibits the same mass gap as seen in the rela- 
tivistic theory and thus the bosonic part of the system still satisfies the Landau criterion for 
superfluidity. Moreover, the phonon dispersion relation retains the same essential form that 
gives rise to the unusual potential for a static probe charge placed in the condensate (2.34). 
Thus the screening of electric charge in the non-relativistic theory is identical to that found 
in the previous chapter. 



3.3 White dwarf cooling 

The presence of charged condensation in the cores of some white dwarf stars could have 
significant observational consequences. In this subsection we follow Ref. [9]. Consider the 
phonon contribution to the specific heat of a white dwarf in which the core has crystallized. 
A conventional phonon in a crystal has a dispersion relation uj oc |k| at low |k|. As a 
result, the contribution of the phonon to the specific heat scales as T^. In contrast, in the 
condensate, the phonon dispersion relation exhibits a mass gap as discussed above. As a 
result, contributions of the bosons to the specific heat at low temperatures arc exponentially 
suppressed as exp(— m^/T). Since for typical densities in helium white dwarfs we have 

~ 3 KeV, the suppression at temperatures below 10^ K will be ~ exp(— 30). 

The electrons exhibit largely the same behavior in both the crystallized and condensate 
phase. At the relevant temperatures they form a degenerate Fermi gas with gapless exci- 
tations near the Fermi surface. As a result their contribution to the specific heat scales 
linearly with temperature. Thus, in the case of the crystallized cores, the contribution of the 
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electrons to the specific heat is subdominant compared to the contribution coming from the 
massless crystal phonon. In the condensate, however, the gapless excitations of the electrons 
will be the dominant contribution to the specific heat. This difference in the specific heat 
between the crystallized phase and the condensate phase has a significant impact on the 
cooling of helium white dwarf stars. 

Using the approach of [34], and following [9] we consider an over-simplified model of a 
reference helium star of mass M = 0.5 Mq with the atmospheric mass fractions of hydrogen, 
and heavy elements (metallicity) respectively equal to 

X ~ 0.99, Zm ^ (0.0002 - 0.002) . (3.12) 

The lower value of the metallicity Z„i ~ 0.0002 is appropriate for the recently discovered 
24 He WDs in NGC 6397 [10], but for completeness, we consider a wider range for this 
parameter. 

It is straightforward to find the following expression for the cooling time of a star in the 
classical regime [34] 
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(3.13) 



where Tf and Tq denote the final and initial core temperatures. The first term in the brackets 
on the right hand side corresponds to cooling due to the classical gas of the ions and the 
second term corresponds to the contribution coming from the Fermi sea. The latter is sub- 
dominant in the range of final temperatures that we are interested in (the factor Z in front 
of this term is due to Z electrons per ion). Since Tf <^ Tq, the age of a dwarf star typically 
doesn't depend on the initial temperature. Neglecting the fermion contribution, we find the 
time that is needed to cool down to the critical temperature Tj = Tc 

where an order of magnitude interval in (3.14) is due to the interval in the envelope metallicity 
composition given in (3.12). We also find the corresponding luminosities 
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L{%) ~ (10« erg/s)^ (^^J ~ 1.5 • (10"^ - 10-^)Lo , (3.15) 

which are in the range of observable luminosities (Lq ~ 3.84 ■ 10"^'^ erg/ s). 

After condensation, the specific heat of the system dramatically drops as the collective 
excitations of the condensed nuclei become massive and "get extinct" . The contribution from 
the Fermi sea, which is strongly suppressed by the value of the Fermi energy, becomes the 
dominant one. The phase transition itself will take some time to complete, and the drop-off 
in specific heat will not be instantaneous. 
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Figure 3.1: A schematic sketch of the luminosity function for hehum-core white dwarfs, 
taken from Ref. [9]. The absolute normalization of the function is set by the constant b 
which is determined by their formation rate. The blue line represents the Mestel regime. 
The shape of the luminosity function near the condensation point depends on the details of 
the corresponding phase transition. 



In the zeroth approximation, we can regard the transition to be very fast, and retain only 
the fermion contribution to the specific heat below Tc. Then, the expression for the age of 
the star for Tf < Tc, reads as follows 
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(3.16) 



Notice the difference of (3.16) from (3.13) - in the former Tf < Tc and it is Tf that enters as 
final temperature in the fermionic part, while Tc should be taken as the final temperature in 
the bosonic part. 

From the ratio of ages, r] = tHe/i'Hei identical helium dwarf stars, with and 

without the interior condensation, we deduce that the charged condensation substantially 
increases the rate of cooling- the age could be twenty times less than it would have been 
without the condensation phase [9]. 

The condensation of the core would induce significant deviations from the classical curve 
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for helium white dwarfs. What is independent of the uncertainties involved in these discus- 
sions, is the fact that the luminosity function (LF) will experience a significant drop-off after 
the charged condensation phase transition is complete. This is due to the "extinction" of 
the bosonic quasiparticles below the phase transition point. In fact, the LF will drop by a 
factor of ~ 200. This is illustrated schematically in Fig. 3.1, taken from Ref. [9]. 

This may be relevant for the results of Ref. [10]. The latter work reported the observed 
termination of a sequence of the 24 He WD's found in a nearby globular cluster NGC 6397. 
The termination takes place before the limiting luminosity is reached. Moreover, it takes 
place before termination is seen in the carbon- and/or oxygen-core white dwarfs observed in 
the same cluster. The direct astrophysical explanation - the disruption of the helium-core 
WD containing binaries by collisions - does not seem to be favored for the local densities in 
the environment where the sequence was observed [10]. It is tempting to speculate that the 
observed termination of the sequence may be a signature of fast cooling due to the charged 
condensation. For these discussions see Ref. [9]. 
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Chapter 4 
Vortex Structure 



4.1 The vortex solution 

In this chapter we argue that the charged condensate has properties somewhat similar to type 
II superconductors. In particular, we show that it can admit solutions that arc similar to the 
Abrikosov vortices [39], originally found in the Ginzburg-Landau model of superconductivity, 
and later recovered in the relativistic abelian Higgs model in [40]. The vortex solution is a 
topologically stable configuration, characterized by a nonzero winding number of the phase 
of the complex scalar field. Asymptotically, the scalar field is given by ~ ve'^, where v 
is the vacuum expectation value of the field and 9 is the azimuthal coordinate. Like the 
Abrikosov vortex, the charged condensate vortex carries a quantized magnetic flux. The 
vortex solution has a higher energy density than the pure condensate solution. However, in 
the presence of a sufficiently high external magnetic field, it becomes energetically favorable 
for the charged condensate to form vortices. 

The obtained vortex line solution exhibits the following structure: it has a narrow cylin- 
drical core where the scalar field changes significantly from a zero to a nonzero value; this 
core is surrounded by a broad halo in which the magnetic fiux is confined. The width of the 
latter region is determined by the penetration depth (i.e., the photon Compton wavelength). 
We refer to the system of the core and the halo as the fiux-tube, or the vortex line. This 
structure is similar to that of the Abrikosov solution. However, unlike the latter, our solu- 
tion also carries a profile of the electrostatic potential within the halo, while this potential 
is exponentially small outside of the halo (i.e., the fiux tube is charge neutral). Hence, in 
terms of differential equations, one has to solve three coupled equations instead of the two 
required in the Abrikosov [39] , or Nielsen-Olesen cases [40] . 

If charged condensation exists in the cores of some white dwarf stars, it is possible that 
vortex structure exists as well. Surface magnetic fields have been found in white dwarf stars 
ranging from 10^ to 10^ Gauss. We will show that the presence of such strong magnetic fields 
should be sufficient for the the creation of vortices in the condensed cores of white dwarfs. 
The existence of vortices in the cores of these stars could have observational consequences. 

The organization of this chapter is as follows: In this section we consider a generic 
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system of charged scalars and oppositely charged fermions in the context of the relativistic 
formahsm presented in chapter 2. We fix the phase of the scalar field to be of the vortex-type 
and solve the corresponding equations of motion. We compare our solutions to those found 
in the abelian Higgs model. In section 4.2 we consider corrections to our solutions due to 
the dynamics of the fermions. In section 4.3 we consider the effects of an external magnetic 
field on the charged condensate and determine the magnitude of the external field for which 
it becomes energetically favorable to form vortices. In section 4.4 we treat specifically the 
case of helium-4 nuclei and electrons. We describe the system in the context of the non- 
relativistic low energy effective field theory developed in the previous chapter. We discuss 
the applicability of the vortex solutions found in this section (4.1) to the helium-4 nuclei 
and electron system. We also consider the effect of a constant rotation on the condensate of 
helium-4 nuclei. 

To find the vortex solution we use the relativistic formahsm developed in chapter 2. We 
assign a charge of +2e to the scalars and — e to the fermions in anticipation of the helium-4 
nuclei and electron system to be discussed later in section 4.4. However, for now we keep our 
considerations general. The scalar field considered below could be any fundamental scalar 
field, or a composite order parameter suited for a problem at hand. Our conclusions are 
independent of the specific charge assignment. The relativistic Lagrangian is given by (see 



1^ = + lid.af + \{2eA^ + /x,5^o - d,af - \ml - eA^J, . (4.1) 



and the photon mass is = 2e ^JJojTjmH. In the following we will consider a system whose 
net charge is zero and thus we will set (Aq) = and /i^ = rnn. 

To find the charged condensate solution (4.4) in the unitary gauge, we fixed the phase 
of the scalar field to zero. We now consider a configuration where the phase is not set to 
zero, nor can it be set to zero everywhere by a non-singular gauge transformation. The 
requirement that the scalar field be single- valued everywhere is satisfied by demanding that 
the change in phase around a closed loop be an integer multiple of 27r. In a system with 
cylindrical symmetry, this is satisfied by setting a = n9, where 9 is the azimuthal coordinate 
and n is an integer. This phase can be removed everywhere by a gauge transformation 
A/^ ^ A/^ + dn{n9), except at the origin where the gauge transformation would be singular. 
The solutions of the equations of motion (4.2), (4.3) where the phase is fixed to a — n9 are 
vortex-type solutions. 



2.8) 



The equations of motion that follow from (4.1) are: 

□ (7 = [(2eA^ + iigSf^o - d^af - mjj] a . 



(4.2) 
(4.3) 



In the unitary gauge the charged condensate solution is given by 




(4.4) 
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At the origin r = 0, in order for the scalar field to be well-defined we must have 

(7 = 0. Here r is the 2D radial coordinate. Far from the origin however, we expect the 
solutions to recover the condensate values (4.4). At large r then, the gauge field takes the 
form 2eAj — > dja, or equivalently Ag — > n/(2er). From this form of the vector potential, it 
follows that this configuration has a quantized magnetic fiux that is related to the integral 
of Aj around a closed loop at infinity: 

^ = y A-dl = y AerdO = ^ . (4.5) 

The magnetic flux is quantized in units of n. The quantization of flux implies the stability 
of the vortex configuration, although it may be possible for a high n vortex to decay into 
multiple vortices of smaller n. 

To solve the equations of motion (4.2), (4.3) for the vortex configuration we switch nota- 
tion to dimensionless variables. The resulting equations are governed by a single parameter k, 
the ratio of the mass of the scalar to the mass of the photon in the condensate: k = mj^/m-y. 
This parameter k is the equivalent to the Ginzburg-Landau parameter in the theory of super- 
conductors which gives the ratio of the penetration depth to the coherence length. For the 
helium white dwarf star, if we take the mass of the helium-4 nuclei to be roughly ttih = 3.7 
GeV and the electron density to be Jo ~ (0.15 — 0.5 Mev)^, then we have k, ~ 10^. In our 
derivations below we frequently take the large k limit. 

We define x = m^r, set Ar = Az = and /ig — ttih, and perform the following change of 
variables: 

myA{x) = 2exAg{x) , (4.6) 
m.,F{x) = 2ea{x), (4.7) 
mHB{x) = Hs + '2eAo{x). (4.8) 
In terms of these new variables equations of motion (4.2), (4.3) become 

d f 1 dA\ 9/ , s / ^ 
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(4.10) 



^^^JJi) -^^B-r. (4.n) 



X dx \ dx J 

The boundary conditions are set by requiring that the solutions asymptote to the condensate 
solutions for large r, while for r = we have Aq = a = dA^/dr = 0: 

Forx^O: A{x)^0, F{x) ^ 0, f ^0. 

Forx^oo: A{x) ^ 1, F{x) ^ 1, B{x)^l. ^ ^ 

We can compare these expressions to those obtained in the usual abelian Higgs model. 
Suppose that instead of Lagrangian (2.8) we had the abehan Higgs Lagrangian: 

= ~Fi^ + kd.^r + k^eA, - d,af - ]{a^ - • (4-13) 
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Using the same change of variables as above and defining — VXv, m^^ — 2ev, the 
equations of motion are: 

^^[-^] = n^-i), (4.14) 



dx \x dx 
I d f dF 

X 



X dx \ dx 



x"^ 



F. (4.15) 



The equation of motion for the vector potential, expressed via A, is the same as in the 
charged condensate model. In the equation for the scalar field, the cr^ term in the abelian 
Higgs model gets replaced in the charged condensate model by a term that depends on the 
electric potential. In addition, in the charged condensate equations the electric potential is 
generally not zero and not constant and has its own equation to satisfy. 

Let us first examine the asymptotic behavior of the solutions to the condensate equations 
for a; — > oo. Far from the origin we expect the fields to be very close to their condensate 
values. Then, on the r.h.s. of equation (4.9), it follows that A(x) — 1 = a{x) is very small. 
If we consider this equation only to first order in small fields then we can approximate the 
scalar field on the r.h.s. of (4.9) as F ~ 1. The solution for A that obeys the appropriate 
boundary conditions is 

A{x) = 1 - CaxKi{x) . (4.16) 

Here Ca is a constant to be determined by the matching of the solutions and Ki{x) is the 
modified Bessel function of the second kind. In the large x hmit this solution for A becomes 

A{x) ^ 1 - Ca^e-\ (4.17) 

To find the asymptotic behavior of B and F we expand these fields in terms of pertur- 
bations above the condensate values, B{x) = 1 + b{x) and F{x) = 1 + f{x), and we assume 
that b{x),f{x) <^ 1. We then substitute these expressions as well as expression (4.17) into 
the equations for B and F and keep only terms linear in the perturbations b{x) and f{x). 
These two equations can be combined to obtain a fourth order differential equation for b{x). 
Using the ansatz b{x) = c^x'^ e^^^ where C6, s, and k are as yet undetermined constants, we 
can find the particular and homogenous solutions for b{x) in the large x limit. We also take 
Ca — 1 which we will justify later. For the particular solution we find 

For the homogenous solution we find s = —1/2 and — {1 ± \J\ — 16/«^)/2. In the limit 
that K is very large, the solution becomes 

bh{x) ^ — [ci sin(v^a;) + C2 cos(\/Ka;)] , (4.19) 

V X 
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for some constants Ci and C2. The complete solution is then 



B{x) = 1 + ^^^2^3^ ^ + [ci sin(v^x) + C2 cos(V^a;)] . (4.20) 

The solution for F{x) can be found once B[x) is known: 

STTTh^ e Ace v^-^ 
= 1 + , 2 I ON [cicos(v^a;) - C2sin(v^a;)] . (4.21) 

Here Ci and C2 are the same integration constants that appear in the expression for B{x). 

As K is large, the second term in the above expressions for B and F dominates the 
asymptotic behavior. For x — > 00 we have 

The asymptotic behavior for the vector potential and the scalar field are similar to that for 
the abelian Higgs model: 

The vector potential A, and thus the magnetic field, are the same in both the charged 
condensate and abelian Higgs models. The asymptotic behavior of the scalar field in the 
abelian Higgs model in the large k limit is dominated by the e~^^ term, as in the charged 
condensate model. 

Notably, this is not the asymptotic behavior for the Abrikosov vortex given in the Niclscn- 
Olsen paper [40]. This discrepancy was first pointed out by L. Pcrivolaropoulos in [41]. The 
incorrect asymptotic behavior is obtained if one similarly expands A{x) as A{x) = 1 + a{x) 
and only keeps terms linear in a{x). This is because the last term in (4.10) and the last term 
in (4.15) are quadratic in a{x) and yet, due to the different exponential dependence of the 
perturbations a{x), b{x) and f{x), these terms can be dominant over terms which are linear 
in 6(a;)and f{x). Linearizing A gives the correct asymptotic behavior of the fields only in 
the limit that k is small. 

The second term in the full expressions for B and F and thus the asymptotic behavior of 
both fields is due strictly to the presence of a nonzero magnetic field. In the absence of any 
magnetic field, the screening of any small perturbation of the fields above their condensate 
values vanishes as e~^^ = Q-^/mm^r ^^^g fQ^^fj 2.25) and (2.34)). 

The above asymptotic expressions are valid as long as x ^ 1, or equivalently r ^ 1/m^. 
We consider now the intermediate region l/y^ ^ a; <C 1, or equivalently 1/M <^ r <^ l/w-y 
where again we have defined M = ^mnm^. At distances much larger than 1/M we assume 
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that the scalar field F is still close to its condensate value. Thus expression (4.16) is still 
valid for A. In this regime then n?' /x^{A — 1)^ ~ -n? /x^. The equations for B and F become 



1 d 

X dx 



dB 
dx 
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X dx \ dx ) x^ 



[1-B^)F. (4.25) 



The solutions are straightforward to find: 



Alternatively, we can once again expand B and F above their condensate values, B{x) = 
1 + b{x) and F{x) — 1 + f{x), and solve for b{x) and f{x). The homogenous solutions for 
h{x) and f{x) are the same as those given above with the same coefficients Ci and C2- Solving 
for the particular solutions gives the full solutions in the linearized approximation: 

B{x) = 1 + + [ci sm{^/^x) + C2 cos{^/kx)] , (4.27) 

lH-^X-^ sjx 

^(^) = 1 - + ~ [''^ cos(x/;^a;) - C2 sin(x/^a;)] . (4.28) 

The coefficients ci and C2 are needed to perform the matching. However, as we'll see below, 
these coefficients will turn out not to be exponentially large, and hence solutions (4.27) and 
(4.28) approximate well the solutions in (4.26). 

The approximations made to find both the homogenous and particular solutions break 
down as x approaches 1/ ^/k. Moreover, f[x) becomes of order 1 at a; ~ ^1 and thus the 
linear approximation in general no longer holds below this scale. 

Finally, we'd like to solve in the r ^ limit. Before we do so, however, wc emphasize 
that validity of this procedure needs some justification. The interparticle separation is given 
by d oc J(7^^^. This corresponds to x oc 1/ At distances shorter than this x we expect 
that an effective field theory would break down and thus it would make little physical sense 
to solve the equations (4.9), (4.10) and (4.11) in this regime. Moreover, the scale 1/M is 
typically shorter than the interparticle separation d, hence, particles at these scales cannot 
in general be modeled by a smooth distribution. 

Both fermions and bosons are in a condensate state in which the location of individual 
particles has uncertainties much greater than the interparticle separation. Hence the latter 
notion loses its meaning as a microscopic characteristic of the system. For this, we'll still 
approximate particle distributions by smooth functions all the way down to the scale ~ 1/M, 
which is a dynamically determined short-distance scale at which weakly coupled expansion 
breaks down [4]. As to solving at scales less that 1/M, we regard this as a purely mathemat- 
ical exercise aimed at finding the matching of the asymptotic solutions for the corresponding 
differential equations for all values of the coordinate x. 
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Taking A, B, and F to be series expansions in small x obeying the appropriate boundary 
conditions, the solutions to (4.9), (4.10) and (4.11) are 

A{x) = aox''-^x\ (4.29) 
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x"^ 



B{x) = bo--, (4.30) 



F{^) = /o 



x-l{K\bl-l) + 2ao) x^ 



(4.31) 



For simplicity we have solved for the case that the winding number n = 1} The coefficients 
oo, /o, and bo, as well as the coefficients Ca, ci, and C2 can be determined by matching the 
above solutions to those in the intermediate region, given by (4.16), (4.27) and (4.28). 

To determine the physically appropriate matching radius, we first use Gauss's law to 
find the charge of the vortex solution. The number density of fermions in the vortex Jq is 
constant and is the same as the number density of fermions in the normal condensate phase. 
We have fixed it so by hand, but will justify this later. The scalar number density is given by 
^JqBF^ and varies as a function of x. Therefore it is not in general equal to its condensate 
value I Jq. The variation of the scalar number density away from its condensate value can 
lead to a net charge density of the vortex core within the vortex halo. In particular, there are 
two competing effects. In the intermediate region I/^/k <^ x ^ 1, both B and F are above 
their condensate values, thus the scalar number density is greater than the scalar number 
density in the condensate. As a; — > 0, however, F ^ and the scalar number density drops 
to zero, significantly below the condensate value. The matching radius should be chosen 
so that these two effects combine to give the appropriate charge density as determined by 
Gauss's law. 

From Gauss's law we can calculate the average charge density of the vortex inside radius 
X = 1. As is usually the case, we can determine the net charge enclosed in a region knowing 
only the form of the potential at the boundary of that region. Equation (4.26) gives the 
potential in the intermediate region independent of matching coefficients Ci and C2. This 
form of the potential, together with Gauss's law, allows us to calculate the net charge of the 
vortex at a; = 1 independent of the matching conditions and the x ^ solutions. 

Gauss's law is given by equation (4.2): 

V^^lo = 2e(2eAo + Us - a)cr^ - e Jo . (4.32) 

The r.h.s. of the above equation is the charge density. Integrating both sides of the above 
expression over the volume of the vortex and dividing by the total volume gives the average 
charge density inside distance x: 

% = 2eJo-^. (4.33) 
V X ax 



^Foi n 1, the leading term in the expansion for F will be oc a;". The expression for B{x) remains 
unchanged and the leading term in the expansion for A{x) = aox^ is also the unchanged. 
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Here V is the volume equal to the length of the vortex times the cross-sectional area and the 
r.h.s. is evaluated at the boundary of the vortex. Using expression (4.26) for B{x) at x = 1, 
the average charge density inside x = 1 is Qcnc/V = —^cJ^/k^. The negative sign indicates 
a dearth of scalars in this region, but, as n is very large, this is a small correction to the 
overall average charge density of scalars oc cJq. To check this result one can hkewise use the 
asymptotic solution for B, expression (4.20), at x = 1. Assuming that the coefficients Ci and 
C2 arc not exponentially large and thus these terms are not dominant in the solution for B 
one finds Qcnc/^ —leJ^jK^. This is consistent with the previous result. Farther 
out, B{x) — \ is exponentially suppressed thus the net charge of the vortex approaches zero 
as X becomes large. 

We can now use this result to determine the matching radius R. Given the smallness 
of the average charge density found inside x = 1, the excess of scalars in the intermediate 
region of the vortex must cancel the shortage of scalars in the x — region to great accuracy. 
Using expressions (4.26) for B and F, it can be shown that this happens when R ~ 
Thus we use this as our matching radius R in what follows. 

We start by matching the solution for A{x) in (4.29) and its first derivative with its 
solution in the intermediate region (4.16). Taking the matching radius to be small, it! <^ 1, 
gives: 

R^ 



7 + ln( :| 



Ca = 1 + ^ , (4.34) 



where 7 is the Euler-Mascheroni constant. As long as R is less than one, is positive. 
Moreover, we see that we were justified in taking ~ 1 in our previous calculations. The 
magnetic field is given by 

/f=?^~. (4.35) 

2e X ax 

Near the origin the magnetic field is of order m^/(2e). For x > it is given by 

w?^Kq{x) / {2e) , where Kq{x) is the modified Bessel function. For x » 1, i.e. for r :§> 
the magnetic field is exponentially small. 

To find the remaining coefficients, we use oq found above and match (4.30) and (4.31) and 
their first derivatives to the appropriate solutions in the intermediate region (4.27), (4.28). 
We now take the matching radius to be = 1/\/k- The solution with the lowest energy is 
one in which the scalar field F(x) is identically zero in the region x < R. The corresponding 
coefficients are 

fco = 1 + ci = «-^/4e(2cos(l)+sin(l)) , 

(4.36) 

/o = 0, C2 - K-^/^e(cos(l) -2sin(l)) . 

In Fig. 4.1 below, the fields are plotted for small r and for k ~ 10^. The radius r = 1/M cor- 
responds to the matching radius x — R — 1/ y/R,. The radius r — d denotes the interparticle 
separation d — Jq . Unlike the magnetic field, the potential and scalar field approach their 
condensate values for x > l/\fK. This is in contrast to the abelian Higgs model in which 
the scalar field is close to its condensate value for x > 1/ i.e. for r > l/mu- 
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Figure 4.1: Small r solutions for the scalar field and electric potential. 

4.2 Fermion dynamics 

In our discussions above we have treated the fermions as a fixed charge background = 
Jo^fio- We relax this assumption now and introduce dynamics for the fermions via the 
Thomas- Fermi approximation. The fermion dynamics are governed by the constant chemical 
potential given in (2.7). The local number density of fermions is determined by the Fermi 
momentum: Jo{x) = p^{x) / {Stt"^) . In this way the number density of the fermions Jq gets 
related to the electric potential Aq. For relativistic fermions 

Mx) = ^ifiF - eAoix)f . (4.37) 

The chemical potential gets fixed by the value of the fermion number density in the con- 
densate phase, where (Ao) = 0. If Jo represents the number density of fermions in the 
condensate, then fip = (Stt^ Jq)^/'^. The photon mass is also defined in terms of Jq: 
= 2e\/ Jo/2mH- In the vortex phase Jo{x) —>■ Jq for large x. 

To include the effects of an x-dependent Jq into our equations, (2.20) gets incorporated 
into the equations of motion (2.9). As a result the equation of motion for B{x) (4.11) gets 
modified. In the linearized equations, the effect is the addition of a new term for b{x) with 
a coefficient which scales as oc mu/^p. However, it turns out that this new term does 
not contribute significantly to the solutions. This is because, in the fourth order differen- 
tial equation for 6(a;), terms with coefficient uih/ixf oc k^/^ are subdominant compared to 
terms with coefficient k"^. Accordingly, the solutions found above in the intermediate and 
asymptotically large regions are still valid. It can be shown that the a; — solutions (4.29), 
(4.30), (4.31) are also unaffected. In physical terms, the inclusion of the fermion dynamics 
via the Thomas- Fermi approximation gives rise to ordinary Debye screening. This screening 
is subdominant compared to other screening effects in the condensate (see [6]). Moreover, 
the profile of Aq[x) away from the core and within the halo is very shallow, giving rise to a 
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very mild dependence of the charge density on x. Hence, the latter can be approximated by 
a constant, as was done in the previous sections. 

As mention in section 2.3, the TF approximation does not capture the possibility of 
exciting gapless modes near the Fermi surface. To include this effect we calculated the 
one-loop correction to the gauge boson propagator, which gives corrections to the static 
potential Aq. We are interested in how this correction compares to the potential found in 
the intermediate region of the vortex (4.27). To estimate its magnitude, we consider a toy 
model of the vortex. We find the potential due to a wire of constant linear charge density Aq 
located at r = 0. The linear charge density of this wire is set by the characteristic charge of 
the vortex: since the scalar charge density varies significantly from its condensate value e Jq 
at scales r < 1/M, it follows that at short distances the linear charge density of the vortex 
can be approximated by 

Ac = 9^. (4.38) 

At large distances the vortex is effectively neutral, as mention above. Thus we expect the 
one-loop contribution to the static potential to be irrelevant at large scales. 
In three dimensions, the charge density of the source is given by 

j^our,,^(r^e,z) = —5ir). (4.39) 
7rr 

The static potential is determined from this source and from the {00} component of the 
gauge boson propagator Doo^ 

Ao(r) = j d^r' Doo{r - r') Jo'°"'"''(r') . (4.40) 

The one-loop correction to the propagator due to this was found in (2.43): 

aem sin(2A;irf) 
^oo(r) = — . (4.41) 
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Here f represents the 3D radius in spherical coordinates, as opposed to the 2D radius r. 

Using this expression together with expression (4.39) in equation (4.40), the correction 
to the static potential is 

, a^Xokl, r sin(2fe^v/i^) 



An upper bound on the potential can be found by taking sm.{2kF\/ z'"^ + r^) 1. After 
integrating, this gives 



^oW < ^^TTTTr^ ^ ~T\ r^- (4-43) 
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On the vortex solution, the leading term in the potential in the intermediate region is 
given by expression (4.27): 

Ao(r) = ^(5(r) - 1) ~ , ^ , . (4.44) 

Given that both kp/mH -C 1 and l/ijUHf) <^ 1, we see that the one-loop correction to the 
potential (4.43) is greatly suppressed compared to the potential found in the vortex solution. 
Thus the excitations of the fermions do not significantly alter the vortex solutions. 

One further effect that we take into consideration is the Landau quantization of the 
fermion energy levels due to the presence of the magnetic field in the interior of the vortex. 
In the presence of an applied magnetic field, the separation between energy levels is given 
by cj = eH/mp where H is the magnetic field and mi? is the fermion mass. Near the core of 
the vortex where H ~ m^/(2e) the separation of levels of the fermions is a; ~ m^/mF- Since 
the photon mass is generally much smaller than the fermion mass, this shift in energy is 
negligible compared to the typical energy of the fermions. 



4.3 Energetics and external fields 

We now turn to the question of when it is energetically favorable to form a vortex in the 
charged condensate. We start by comparing the average energy density of the vortex to the 
energy density of the pure condensate. Above, using Gauss's law, we found that inside the 
distance x — 1 the vortex has a small negative charge density implying that in this region 
the average scalar number density is lower than in the condensate phase. At distances a; 3> 1 
this charge density is exponentially suppressed indicating that the net charge of the vortex 
is zero and thus the total average scalar number density is the same in both the vortex phase 
and the condensate phase. In calculating the average energy density of the vortex inside the 
distance x = 1, we are not interested in the contribution to the energy due to the discrepancy 
in the number of scalars between the vortex phase and the condensate. This contribution 
to the overall difference in energy vanishes at large distances. Thus we calculate the energy 
density of the system using Ti' = Ti. — figJ^^'^^'^^- The additional term effectively subtracts off 
the energy density due to the scalar number density. We compare Ti' in the vortex phase to 
H' in the condensate. 

The Hamiltonian density H' can be calculated from the Lagrangian C (4.1): 

n' = + ^E^ + l{2eAo + i^s - a)V' - //,(2eAo + - q;)(t' , (4.45) 

We have simplified the Hamiltonian using the equations of motion (4.2), (4.3) and have taken 
boundary terms to be negligible. The magnetic field H and the electric field E are defined 
as usual 

H^^-^{rAe), E ^ . (4.46) 
r dr ar 

The fourth term in the Hamiltonian is exactly —^gj^'^'^^'^^ as we would expect. The third 
term is due to the energy of the scalar field. Unlike in the abelian Higgs model, the energy 
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density of the scalar field in the center of the vortex, i.e. in the "normal" phase, is lower than 
in the condensate phase. However, the energy density of the vortex is still greater than that 
of the condensate alone, due to the gradients of the scalar field and due to the intermediate 
region 1 / s/k -C a; -C 1 in which the values of both the potential and the scalar field are 
greater than their condensate values. This contribution to the energy density is roughly 
equal in magnitude to the contributions coming from the electric and magnetic fields. On 
the condensate solution, the Hamiltonian density is identically zero: Ti'cc ~ 0- 

For large deviations away from the condensate arc exponentially suppressed and thus 
differences in energy between the two phases are negligible. So to find the average energy 
density within the vortex, we integrate the Hamiltonian density over an area of radius x — \ 
and then divide by the total area. The average energy density within the radius x < 1 is 

eave= I IxdxV! . (4.47) 







The Hamiltonian density can be further simplified using the equations of motion (see Ref . [42] 
for more details). In terms of the dimensionless variables defined above equation (4.47) 
becomes: 

1 

Cave = ^mnJo / X dx {B{x) - 1) (3 - F{xf) . (4.48) 
^ Jo 

For the region x < we use solutions (4.30) and (4.31) and in the intermediate region 

l/yK < a; < 1 we use solutions (4.27) and (4.28) with the coefficients found from matching 
(4.36). Upon integration, the average energy density is 

eave = ^(log«:+14) . (4.49) 

The numerical coefficients should not be taken too literally given the approximations made 
in obtaining the solutions which yield the above result. However, the overall scaling of the 
energy density eave oc mj^ Jo(log ti)/ is remarkably independent of the matching radius and 
other details of the solutions. As the energy density of the condensate is effectively zero 
{V!qq — 0), the above expression represents the difference in energy between the two phases. 

To see when it is energetically favorable for the condensate to form vortices, we now 
consider placing the condensate in an external field i^cxt pointed along the 2;-axis. We shall see 
that the magnetic properties of the charged condensate resemble those of a superconductor. 
In particular, when k ^ 1, the charged condensate resembles a type II superconductor. 
When an external magnetic field i?ext is apphed to the condensate, below a critical value 
Hci surface currents oppose the penetration of the field and the induction is zero in the 
bulk of the condensate. For iJcxt > H^i magnetic flux penetrates the condensate in the form 
of vortices. At another critical value of the magnetic field the normal phase is restored 
and the induction B^^^ is equal to the applied field i^ext- In what follows we determine the 
critical values of the fields Hd and Hc2- 

Given the energy density e of the vortex phase above the pure condensate phase, we 
can find the value of the magnetic field Hd at which it becomes energetically favorable 
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to form vortices. In the absence of an external field, it is never energetically favorable to 
form vortices as the energy density of a vortex is greater than that of the pure condensate. 
In the presence of a small external magnetic field, below Hd, the condensate must expel 
the magnetic field entirely from its bulk in order to remain in the condensate phase. This 
requires energy; the energy per volume needed to expel the external field is l-f^ext- vortices 
form in the condensate then the energy required to expel the magnetic field is smaller than 
if the field were to be completely expelled. More specifically, if the vortices give rise to an 
average magnetic field in the condensate -Bind, then the energy needed to expel the remaining 
magnetic field would be |(i?ext ~ -Bind)^- The energy gained by forming vortices is the 
difference between this energy and the energy required to expel the magnetic field entirely. 
Assuming that Sind is small compared to i?ext near the transition point, this difference can 
be approximated by ^(2Bi^^Hcxt) ■ Thus for a high enough external field, the energy e lost in 
creating a vortex is compensated by the energy gained in expelling a smaller magnetic field 
BindHcxt- In order for formation to be energetically possible, we must have e < BindH^xt- 
The equahty determines the critical external field Hd- (See Ref. [39].) 

Suppose the number of vortices per area in the condensate is given by N. Then the 
energy density due to the formation of vortices is given by e = NX where A is the energy 
density per unit length of a vortex. Using eave found above (4.48) as the energy density of a 
single vortex. 



A=^eave. (4.50) 
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The induction Bi^id is given by 

And-A^^A-dl=^. (4.51) 
Combining these expressions, the critical field Hd = A^A/(i?ind) is given by 

2 

tn 

Hci = ^ (log(«;) + 14) . (4.52) 
8e 

The final expression for Hd is independent of the number density of vortices N. It follows 
that if it is energetically favorable to create one vortex, then it will be even more energet- 
ically favorable to create many, up to the point than interactions between vortices become 
significant. At distances greater than r — we expect fields outside the vortices to be 

exponentially suppressed and thus the vortices to be effectively non-interacting. So at the 
transition point Hd, it is likely that the number density of vortices is of the order ~ m^/TT. 

If we take Jo — (0.15 — 0.5 MeV)^, a reasonable value for white dwarfs, this gives a 
magnetic field of roughly Hd — (10^ — 10^) Gauss. Thus, the vortex lines should be expected 
to be present in the bulk of the helium-core white dwarf stars with strong enough magnetic 
fields. The condensed-core WD's with fields of strength less than ~ 10^ Gauss would expel 
their magnetic fields. Such stars are expected to constitute a significant fraction of the 
helium-core WD's. 
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A sufficiently high magnetic field will disrupt the condensate entirely. One way to ap- 
proximate the magnetic field at which this transition occurs is to consider the density of 
vortices in high external magnetic fields. When the cores of the vortices begin to overlap, 
then the scalars are mostly returned to the normal phase. We define the core of the vortex 
to correspond to x = or equivalently, r = 1/M as this is the region in which the VEV 

of the scalar field drops to zero. As N is the number of vortices per area, at the transition 
point N ~ M^/tt. We define the critical external field at this point to be Hc2- When the 
condensate enters the normal phase, the induction B^^^ is equal to the external magnetic 
field Hc2 = -Bind = 27rA^/(2e). From these two expressions we find Hc2- 

i/c2 = — . (4.53) 
e 

For Jo ~ (0.15 - 0.5 MeV)^, Hc2 ^ (10^^ - 10^^) Gauss. This is well above the values of the 
fields expected to be present in a majority of white dwarf stars. Thus the external magnetic 
field is unlikely to be large enough to push the condensate into the normal phase. 

It should be noted that both Hd and Hc2 given above were determined at zero temper- 
ature. Generally, we expect these expressions (4.52), (4.53) to be vahd at temperatures well 
below the condensation temperature. 

Finally, in type II superconductors the dependence of the critical temperature on the 
magnetic field is well-approximated by T'^ /T"^ ~ {He — H)/Hc, where is the transition 
temperature when the magnetic field H is present. We expect a similar relation to be valid in 
our case too. Hence, as long as the value of the magnetic field is not too close to either critical 
value, the change of the transition temperature due to the magnetic field should be small. 
Near the critical values, however, the change of the phase transition temperatures (from the 
normal to the vortex phase and from the vortex phase to the phase with no magnetic field) 
could change significantly. The would be crystallization temperature will also change, and 
the charged condensation may or may not be favorable for close-to-critical magnetic fields. 



4.4 Comments on white dwarf stars 

To appropriately describe vortices in the cores of helium white dwarf stars, we should use 
the low energy effective Lagrangian description of the charged condensate developed in the 
previous chapter: 

A// = V (^{2eAo - doT)^' - ^(V.S)^ - ^(2eA, - d.TfE^^ . (4.54) 

In what follows we consider how this alternative formalism changes the vortex solutions 
found above. 

It was found in section 3.2 that, to linear order, the equations of motion for small per- 
turbations above the condensate were the same in the low energy effective theory as in the 
relativistic effective theory, up to a time derivative for the scalar perturbation r. As we 
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are concerned with static solutions, this time derivative will not change the vortex solu- 
tions in the low energy theory. In order to find vortex solutions in the intermediate region 
<^ X <^ 1 and in the asymptotic region a; ^ 1, we treated the electric potential 
Aq and the scalar field a in the linear approximation, but kept higher order terms for the 
vector potential. Thus to determine the applicability of the solutions found above to the 
non-relativistic effective theory, we should consider higher order terms in Aj than the ones 
given in (3.9). We also restore the phase F. The equations of motion to next-to-leading-order 
that follow from (4.54) are 

-9^F^o = 2e [1 + 2C2E2(2eAo - doT)] - eJo , (4.55) 
-d^Fi^j = 2e(2eAj - djr)T,^ , (4.56) 
-V^E = [2mjy(2eA - doT) + AC2mH^'^(2eAo - doT)^] E - (2e^j)^E . (4.57) 

If we take S = ^nino and C2 = l/(m//Jo), then the first two equations of motion above 
(4.55), (4.56) are the same as in the non-relativistic case (4.2), up to second order in small 
fields. The third equation (4.57) has an extra factor of (2e74o)^E compared to equation 
(4.3). However, since this term is second order in Aq and we treated Aq in the linear 
approximation, this does not alter our solutions in the intermediate and asymptotic regions. 
Thus for C2 = l/(m//Jo), the vortex solutions found above for x ^ are also solutions 

for the non-relativistic effective theory. 

As in the relativistic case, the solutions formally break down near x — l/y/n when 
the change in the scalar field becomes of order 1 and thus the linear approximation is no 
longer valid. More realistically though, we do not expect the effective field theory to hold 
at distances shorter than the interparticle separation x oc The effective field theory 

will cease to be a valid description of the physics before reaching x = l/y/R. Moreover, we 
do not in general expect that the low energy effective theory will obey the Lorentz invariant 
condition C2 = l/(mijJo). Instead, C2 must be fixed by the particular physics of the system. 

We can use the non-relativistic formalism to consider the effects of the rotation of a white 
dwarf star on the helium-4 nuclei and electron system. In this formalism the scalar number 
density and current density are given respectively by: 

jscalax ^ ^ jscalar ^ [{Dj^)*^ - ^*{Dj^)] . (4.58) 

The number density is related to the current density by J^'^'^^'^^ = J^'^'^^'^'^Vj where vj is the 
velocity vector of the rotating scalar particles. Using the change of variables defined above, 
$ = Eexp(ir), we can use these expressions to find Vj: 

V. = — (2eA. - ar) . (4.59) 
ruH ^ ■< ' 

This known result is notably different from that of a superfluid in which the scalar field does 
not couple to a gauge field. In the absence of the Aj term in the above expression, one would 
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conclude that V x v = and thus the scalar condensate does not support rotation. Instead, 
in the presence of the gauge field we find 



Vxv=— H. (4.60) 
rriH 

The magnetic field H is called the London field [43]. 

The velocity vector v can be written in term of the angular velocity v = f2 x r. It follows 
that, for constant ft, the rotation of v is given by V x v = 217. Accordingly, the magnetic 
field can be expressed in terms of the angular velocity: 



H = ^f2 = ^n. (4.61) 
2e m2 ^ ^ 



7 



Here Jo is the fermionic number density. Thus the condensate of helium-4 nuclei can rotate 
with the rest of the star, unlike a neutral condensate. The consequence is a small, constant 
magnetic field in the bulk of the condensate. 

Varying the Lagrangians (3.5) and (3.6) with respect to Aj gives 2eJj"^^^^'^' = eJj, where 
Jj is the fermion current density. Using Jj = JoVj, it follows that the fermion velocity vector 
is equal to the scalar velocity vector. The electrons and the helium-4 nuclei rotate together 
in the core of the star. At the surface however, there is a thin layer of hehum-4 nuclei that 
is slightly out of rotation with the rest of the star. This feature becomes evident upon finite 
volume regularization of the system. The thickness of the layer is roughly This surface 

layer is what gives rise to the London field in the interior of the star [44]. To estimate the 
angular velocity of a helium white dwarf star we take Q ~ 10^^ Hz. The resulting London 
field is if ~ 10~^ Gauss. This field is present even in the absence of vortices. However, it is 
too small to affect any of the results given above. 

The existence of vortices in white dwarf stars with condensed cores could have observa- 
tional consequences. In particular, one could consider the scattering of light near the surface 
of the star. Moreover, the discussions in the previous chapter concerning the specific heat 
of the star would be modified with the existence of vortices, and calculations of the cooling 
rate and the luminosity function would have to be adjusted accordingly. 
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